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ATTEMPTS TO DETECT RESONANCE SCATTERING IN Zn™; THE EFFECT 
OF ZERO- POINT VIBRATIONS* 


R. V. Pound and G. A. Rebka, Jr. 
Lyman Laboratory of Physics, Harvard University, Cambridge, Massachusetts 
(Received March 23, 1960) 


The 93-kev y ray from the 9.4-ysec level of 
Zn” would appear to offer an extreme example 
of the fractionally narrow y rays in solids dis- 
covered by Mossbauer.’ It has been suggested 
as an ideal medium for studying such small fre- 
quency shifts as would result from, for example, 
gravity.2*>* We wish to describe some unsuccess- 
ful attempts to detect resonant absorption of this 
y ray and to offer an explanation which has more 
general applicability and complicates the exten- 
sion of the technique much beyond the narrowness 
presently achievable with Fe*”’.*~” 

The 93-kev y ray is emitted as the final step in 
the decay of 78-hr Ga®™ to stable Zn®™. Ga” was 
prepared with good radiochemical purity by (a, 2n) 
reaction in Cu®® by about 10 microampere hours 
of bombardment by the MIT cyclotron. 

Resonant absorption experiments were attempted 
for one combination of source and absorber at the 
temperature of liquid nitrogen and, for this and 
three other combinations, at the temperature of 
liquid helium. The experiments at helium tem- 
perature used a cryostat with a thin bottom to 
allow the 93-kev y ray to be measured externally. 
In all experiments comparison was made of the 
y-ray intensity transmitted through the absorber 
with the source stationary to that with it vibrating 
at an rms velocity of about 10-° cm/sec, about 
one hundred times larger than the expected line- 
width. The vibration was produced by a quartz 
electromechanical transducer which was also 
tooled. In all cases a null effect was found to the 
accuracy of the statistical uncertainty of + 0.1% 


of the counts in one channel. Correcting for back- 
ground with the tenfold nonresonant absorption 
normally present in the absorbers leaves the re- 
sult that no resonant absorption greater than 

0.3% was present in any of the combinations. 

The combinations tried were: 

1. Source Ga® dissolved in 95% aluminum 5% 
by weight zinc alloy, and absorber 90% aluminum 
10% zinc by weight, 4 cm thick. It was hoped 
that the aluminum would provide a cubic environ- 
ment to avoid electric quadrupole splitting as 
well as a moderately high effective Debye tem- 
perature. Unfortunately the necessarily high zinc 
concentration probably destroys lattice symmetry 
seriously enough to broaden the lines through the 
inhomogeneous electric quadrupole interaction. 
The 11-kc/sec linewidth in the absorber result- 
ing from the 9.4- usec half-life is comparable to 
widths obtained in nuclear paramagnetic reso- 
nance and it is known that even 0.5% of zinc 
broadens the Al magnetic resonance line.* One 
would not expect the Zn to be less affected. Ig- 
noring that broadening and supposing, optimisti- 
cally, that the 375°K Debye temperature of the 
Al host applied to the recoil dynamics of the zinc, 
one would expect an effect of 0.4% at nitrogen 
temperature and 2.2% at helium temperature. 

2. Source Ga dissolved in polycrystalline 
zinc and absorber 5-mm thickness of polycrys- 
talline zinc. In this case one would expect a re- 
duction of the cross section by a factor of 13/50 
because of electric quadrupole splitting in the 
hexagonal close-packed lattice. The Debye tem- 
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perature of zinc is strongly temperature depend- 
ent. If one uses a value of 240°K, one expects 
only 0.08% effect, but if one uses the value of 
305°K characteristic of the specific heat at 4°K 
or of the elastic constants, one expects 0.56%. 
Undoubtedly this estimate is overly optimistic 
because the low-temperature specific heat meas- 
urements and the elastic constants relate only to 
the low-frequency end of the spectrum of lattice 
waves whereas the high-frequency part dominates 
the recoil dynamics. 

3. Source Ga™ dissolved into a single crystal 
of zinc with a 5-mm-thick single zinc crystal 
absorber. These crystals were oriented with c 
axes perpendicular to the direction of propaga- 
tion of the y rays observed and parallel to one 
another. It was hoped that the relative hardness 
of zinc in the direction perpendicular to the c 
axis would enhance the recoil-free component of 
the y ray above that predicted from the Debye 
temperature.’ Again, this is mainly connected 
with the low-frequency part of the phonon spec- 
trum. In addition, the quadrupole splitting would 
reduce the cross section less with oriented single 
crystals than with polycrystalline samples. 

4. Source Ga left in the cubic copper foil but 
annealed to relieve radiation damage, and about 
one-centimeter thickness of natural sphalerite, 
ZnS, from Franklin, New Jersey, reputedly of 
relatively high purity, used as absorber. In this 
case broadening and reduction of cross section 
by electric quadrupole coupling should have been 
avoided if the sphalerite was adequately pure. 
The moderately high Debye temperature of Cu 
should prove helpful. Using 315°K for the Debye 
temperature of copper and 270°K for ZnS leads 
to an expected effect of 1.15%. 

The failure of the above attempts may involve 
several causes. The effects of the electron cap- 
ture that converts Ga® to Zn™ and of the y rays 
leading to the 9.4- sec state may not be healed 
adequately quickly and completely at helium tem- 
peratures to avoid broadening by electric quad- 
rupole interactions. More fundamental, however, 
is an effect related to that involved in the pre- 
dicted’®, *! and observed’ variation with tempera- 
ture of the frequency of the y ray of Fe”. The 
fractional width of the unbroadened line in Zn™ 
should be only 10-**. It is now clear that the 
lattice energy shifts the frequency of the y ray 
through the second order Doppler effect and that 
when different lattices are used for source and 
absorber the y-ray and absorption line will be at 
different frequencies. Large shifts compared to 
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10-*5 have been found between source and ab- 
sorbers of Fe*’ even when both are at a common 
temperature and ostensibly the same material,” 
Differences of Debye temperature were suggeste 
as the cause. Differences of Debye temperature 
should result in such relative shifts even at the 
absolute zero of temperature owing to the differ- 
ence in zero-point energies. Assuming a Debye 
spectrum for the lattice oscillators, the zero- 
point energy of 9k6)/8 per atom leads for Zn” 
to a fractional decrease of frequency of 0.772 
x10-** per degree increase in Debye temperature 
If the zero-point vibration of the Zn™ in copper 
has an effective Debye temperature 45°K higher 
than that of the Zn™ in ZnS, then the absorption 
line in the ZnS would be fractionally 35 full line- 
widths higher in frequency. 

Resonant absorption could be observed only by | | 
compensating that shift by a constant relative 1 
velocity. To find the line without prior knowledg| , 
of the displacement would require an extensive 
search and a good velocity spectrometer. Fur- 1 
thermore, inhomogeneity of the lattice could in- ‘ 
troduce important broadening and consequent 
weakening of the lines through small variations 
of the zero-point energies. Only in the experi- 
ments using zinc directly, of those reported 
above, would there appear to have been hope of 
finding absorption by the method used. For those 
cases the effect expected was probably too small 


a. = ase 25 oe 




















to be observed. I 
Addition of the effect of the zero-point vibra- ( 
tions to the shifts induced thermally modifies the} 
values of the differences in 6p needed to explain — 4 
the shifts at a common temperature between ab- 3 
sorber and source reported previously for Fe"." Fo 
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SEARCH FOR THE ANISOTROPY OF INERTIA USING THE MOSSBAUER EFFECT IN Fe’? 


C. W. Sherwin, H. Frauenfelder, E. L. Garwin, E. Liischer, S. Margulies, and R. N. Peacock 
University of Illinois, Urbana, Illinois 
(Received March 18, 1960) 


The width of the central resonance absorption 
line in the Méssbauer effect in Fe®’ has been 
found by several observers to be considerably 
greater than the natural linewidth.’ In all 
these cases the emitting and the absorbing nuclei 
were randomly oriented, and the possibility 
arises that this widening is due to an external 
perturbation which affects the nuclear levels. 

An intriguing mechanism for line broadening 
has been pointed out by Cocconi and Salpeter.*»® 
The locally asymmetric distribution of matter 
inthe universe, in accordance with Mach’s 
principle, may cause a local anisotropy 4M 
=4(M’-M”) in inertia. Here M’ is the inertial 
mass for matter accelerated towards the center 
of the Galaxy and M” the inertial mass for ac- 
celeration perpendicular to this direction. Such 
an anisotropy AM would give rise to a shift in 
atomic’ and nuclear® energy levels. In the case 
of Fe’, the shift would occur only for the ex- 
cited state with spin 3/2. Cocconi and Salpeter 
compute the level shift AE to be (AM/M)TP.,, 
where T is the average kinetic energy of the 
nucleon responsible for the transition, and P, 
has the value 1/5 for the most favorable relative 
orientation between the nuclear angular momen- 
tum and the line to the galactic center. Assum- 
ing T=10 Mev and a line broadening of 10° ev, 
Cocconi and Salpeter conclude that AM/M has 
an upper limit of about 1 part in 10". 

Since the excess broadening of the center reso- 
tance line is large, it is easy to determine 
whether it arises from effects due to the Galaxy. 
hour first test, we placed both the Fe*”’ source 
and the Fe absorber in parallel magnetic fields 
of about 1000 gauss. In this condition, the atomic 
magnetic fields at the emitting and absorbing 
welei are parallel and the postulated galactic 


shifts should be the same for each of the cor- 
responding magnetic sublevels in source and ab- 
sorber. Thus any line broadening due to ani- 
sotropic inertia should disappear. The experi- 
ment showed, however, that over a twenty-four 
hour period there was no significant change in 
the linewidth between the two conditions: (1) no 
magnetic field, random domain alignments, and 
(2) source and absorber in parallel magnetic 
fields. The accuracy of this test is limited by 
the fact that, although the domains are aligned, 
the local magnetic fields at the nuclei in the 
source and absorber are different, producing 
an additional broadening effect. The results 
showed clearly that the majority of the excess 
broadening was not due to the Galaxy. 

To perform a more sensitive test for a possi- 
ble galactic effect, a magnetic field of about 600 
gauss at the moving horizontal absorber*® was 
aligned in the north-south direction. The source 
was placed in a magnetic field of similar mag- 
nitude, but lying in a vertical, north-south plane, 
and oriented at +45° with respect to the horizon- 
tal (see Fig. 1). Each sidereal day, the galactic 
center rotates once around a line parallel to the 
earth’s axis, changing its direction with respect 
to the magnetic fields of source and absorber. 
Cocconi and Salpeter assume that the level shift 
reaches an extreme value when the magnetic 
field H is perpendicular to the line toward the 
galactic center, and again when it is parallel to 
that direction. At an intermediate angle the level 
shift vanishes. 

In the experimental arrangement, shown sche- 
matically in Fig. 1, there is an interval of about 
10 hours, centered at 17:30 sidereal time, when 
the line toward the galactic center makes nearly 
equal angles with the two magnetic field direc - 
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tions. Under these conditions, no relative shift 
is expected in the upper state energy levels be- 
tween the source and the absorber. In contrast, 
there is about a 3-hour interval, centered at 
06:30 sidereal time, when the galactic line is 
nearly 90° with respect to the absorber, but 
nearly 45° with respect to the source. During 
this latter period, the source should be essen- 
tially unaffected by the Galaxy, but the energy 
levels of the upper state in the absorber should 
experience their maximum shift. With the Gal- 
axy at the nadir, the shape of the center line 
should change, since it splits into two lines near- 
ly symmetrically displaced in frequency, but 
whose intensities have a ratio of about 2:1. 

We have made a number of runs through a 
complete sidereal cycle using the arrangement 
of Fig. 1 looking for periodic variation in (a) the 
amplitude of the center of the line, and (b) the 
amplitude of four Doppler-shifted points on the 
side of the line. We also examined the shape of 
the entire center line at various sidereal times 
to see if there were any obvious deviations from 
symmetry. 

In addition, we made observations with the 
angle ¢ of Fig. 1 at -45°. For this arrangement, 
the line to the Galaxy makes equal angles with 
the two magnetic fields twice a day (at about 
13:00 and 22:00, sidereal time), when no rela- 
tive shift between the excited levels of the source 
and the absorber is expected, but a near-maxi- 
mum relative shift is expected when the Galaxy 
is near the nadir (05:00) and also when it is due 
south (17:30). 

The experiments were performed during the 
period March 2 through March 12, 1960, at lat. 
40°N, long. 88°W. 

Using the arrangement of Fig. 1 for ¢ = 45° 
and for ¢ = -45° the amplitude of the center point 
of the absorption line was constant as a function 
of sidereal time to within the expected statistical 
fluctuations of +0.3% with the confidence level 
of 40% by the x”? test. Similarly, the amplitude 
at the half-intensity point was constant to within 
+0.3% with a confidence level of 25%. 

As an additional check, we looked for fluctua - 
tions in the amplitude of the center point of the 
line for the case where the magnetic fields were 
parallel, both horizontal and pointing north, 
@=0. For these conditions no variation is ex- 
pected and no significant time variation was 
found. 

With the experimental linewidth of 1.2 x10 
ev, a splitting of the center line of magnitude 
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FIG. 1. The geometrical relationships between the 
magnetic fields on the source and the absorber with 
respect to the direction to the center of the Galaxy. 
Tests were also made with $=-45°. 


d 
AE =1x10~ ev would have resulted in a clearly 
observable sidereal variation. Thus, since 1 
E =14.4 kev, AE/E is less than 1x10™. 
In order to calculate AM/Mfrom AE/E certain ( 
assumptions must be made. Following Cocconi \ 
and Salpeter, and taking T = 10 Mev and P,=1/5, | { 
we get the limit AM/M<5x107'*. Both assump- [ | 
tions are, however, open to debate, especially 
since the magnetic moment of the ground state’ I 
and of the excited state® of Fe®’ are considerably [ i 
smaller than one would expect from the shell f 
model. d 
a 





Supported in part by the joint program of the Office c 
of Naval Research, the U. S. Atomic Energy Commis- te 
sion, and by the U. S. Air Force through the Air Force — f 
Office of Scientific Research and through the tri-Ser- i 
vice contract with the Coordinated Sciences Laboratory 
administered by the U. S. Army Signal Corps. 

'R. V. Pound and G. A. Rebka, Jr., Phys. Rev. 


Letters 3, 554 (1959). p 
2S. S. Hanna, J. Heberle, C. Littlejohn, G. J. Per- 
low, R. S. Preston, and D. H. Vincent, Phys. Rev. q 


Letters 4, 28 (1960). 

3g. DePasquali, H. Frauenfelder, S. Margulies, 
and R. N. Peacock, Phys. Rev. Letters 4, 71 (1960). 

‘T. E. Cranshaw, J. P. Schiffer, and A. B. White- 
head, Phys. Rev. Letters 4, 163 (1960). 

5G. Cocconi and E. E. Salpeter, Nuovo cimento 10, 
646 (1958). 

6G. Cocconi and E. E. Salpeter, Phys. Rev. Letters 
4, 176 (1960). 

'G. W. Ludwig and H. H. Woodbury, Phys. Rev. 


















of 
17:30 
ime 


| the 
ith 


arly 


rtain 
‘oni 
1/5, 
imp- 
lly 
ate 
rably 


)ffice 
nmis- 
Force 
Ser- 
ratory 


hite- 





VYotuME 4, NUMBER 8 





PHYSICAL REVIEW LETTERS 








APRIL 15, 1960 








117, 1286 (1960). 
“ts, §. Hanna, J. Heberle, C. Littlejohn, G. J. Per- 






low, R. S. Preston, and D. H. Vincent, Phys. Rev. 
Letters 4, 177 (1960). 





HIGH -FREQUENCY STUDIES ON SUPERCONDUCTING TIN* 


M. S. Dresselhaus! and G. Dresselhaus 
Laboratory of Atomic and Solid-State Physics and Department of Physics, Cornell University, Ithaca, New York 
(Received March 17, 1960) 


Microwave studies of normal metals both with 
and without a static magnetic field have provided 
information on the shape of the Fermi surface’ 
and on the effective masses of the electrons.” In 
this note, the surface impedance measurements 
at 1 kMc/sec on superconducting tin*® both with 
and without a magnetic field are analyzed to study 
changes in the band structure which accompany 
the transition from the normal to superconducting 
metal. It is found that for superconductors the 
zero-field experiments are relatively insensitive 
tothe band structure as compared with experi- 
ments on the field variation of the surface im- 
pedance. Secondly, it is suggested that different 
carriers are important in the two problems, 
with the “heavy” electrons dominating the zero- 
field measurements, and the “light” electrons 
more important in the field-dependent studies. 

The recent extension‘ of the Serber two-fluid 
model calculation of the superconducting surface 
impedance Z to include the term proportional to 
H,’ in Z(H,) allows an analysis to be made of the 
data for H,>0, H, being the static magnetic field 
atthe surface. The temperature dependence of 
the normal and superconducting electron con- 
centrations in the superconducting state is taken 
to be that predicted by the Gorter-Casimir two- 
fluid model, i.e., N,, =Nof* and N,=No(1 -¢*), 
in which No, N,, and N, are the total, normal, 
and superconducting electron concentrations, 
respectively, and ¢ = T/T. is the reduced tem- 
perature. 

For both Hg =0 and 0<Ho<H,, H, being the 
critical field at the temperature T, the Serber 
two-fluid treatment predicts the magnitude and 
‘emperature dependence of Z in terms of the 
parameters (m*y F) and A,, which characterize 
agiven specimen. The quantities m* and vr 
are the effective mass and Fermi velocity of the 
tormal electrons in the superconducting state 
and 4. is the penetration depth A of the dc mag- 
letic field evaluated at T=0°K. These quantities 
determine a dimensionless parameter a, given 
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in which 6, is the superconducting skin depth for 
the rf fields in the extreme anomalous limit ap- 
propriate to the total electron concentration N,. 
The quantities N, and A, are known from other 
measurements. The temperature dependence of 
a, is attributed to the variation of (m*v,) with T. 
The value of a, =0.16 is appropriate for Sn, if 
m* =m, (free electron mass), vp = 6.610" cm/ 
sec (normal state value)° and A, = 5.01075 cm.® 
The determination of a, for a single-crystal tin 
sample* was made for both H,=0 and H,>0. For 
comparison, the field variation of Z was also 
measured and analyzed for the same sample in 
the normal state.” Figure 1 shows a reduced 
temperature plot of the zero-field surface re- 
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FIG, 1. Plot showing zero-field surface reactance 


vs reduced temperature. The solid curves give the 
dimensionless quantity % = |c?X(0)/16vA] as calculated 
from the two-fluid model for three values of ay (see 
text). The experimental points are for a single- 
crystal tin sample at 1 kMc/sec. 
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actance X(0) in dimensionless units, % =[c?X(0)/ 
16mvaA], as calculated from Serber’s model for 
a@,=0.16, 0.28, and 2.17. The quantities v and 
c are the frequency and velocity of light, re- 
spectively. It is seen that the experimental 
points best fit a,=0.16. For comparison, sur- 
face impedance measurements on this sample 
in the normal state yield an effective mass 
(m*/m,) = 0.6 and a relaxation time 7 =9 x10“ 
sec. The curve a, =0.28 in Fig. 1 was deter- 
mined from the measurements of the frequency 
dependence of the surface resistance of Sn by 
Kaplan, Nethercot, and Boorse® in the frequency 
range between 17 and 77 kMc/sec. The disagree- 
ment between the curve for a, = 0.28 and the 
experimental points indicates that the two-fluid 
model analysis used here is not applicable to the 
mm frequency range. However, at frequencies 
v<1.5 kMc/sec a two-fluid model is expected to 
be adequate for Sn.° 

Figure 2 shows a theoretical plot as a function 
of reduced temperature of the longitudinal field 
variation of the surface reactance AX 7 (HQ) 
=[X 7 (Hg) -X(0)] for a, =2.17, which fits the ex- 
perimental points in the region ¢<0.7. At higher 





AX, in 10°n 
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FIG, 2. Plot of field variation of surface reactance 
vs reduced temperature. The solid curve is calculated 
for @)=2.17. The experimental points give the term 
in H,? of AX (A) for a tin sample at 1 kMc/sec, t<0.7. 
AXL corresponding to a)=0.16 is too small by a factor 
107". The value of Hy used to evaluate AX 7(H)) is the 
critical field at the temperature ¢. 
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temperatures, ¢>0.7, the data are fitted by a, -5 
In terms of the parameter y = (m*vp/movp’), 
up’ being the Fermi velocity in the normal state, 
a@,=2.17 and a,=5 correspond to y =3.8 x10 
and y =3x1075, respectively. 

In the presence of H,, both the magnitude and 
the temperature dependence of AX ; (Hg) are de- 
termined by (m* vp). The magnitude of AX 7 (Ho) 
is particularly sensitive to the choice of a@,, with 
AX ; (Hg) <a@.° On the other hand, the zero- 
field experiments are relatively insensitive, a, 
being determined by the temperature at which 
% departs from the constant value of 7/2. To 
emphasize the fact that different values of a, are | | 
needed to fit the data for H, =0 as opposed to 


H,>0, a curve of x; is plotted in Fig. 1 for a, ( 
=2.17. Furthermore, the parameter a, =0.16 
which fits the zero-field data yields a value of : 


AX, too small by a factor of 10’. 

Anisotropy in the band structure is indicated 
by a comparison between the longitudinal and 
transverse field experiments. The transverse 
field data for X are fitted by a, =2.29 in the low- 
temperature range and a, =6 in the temperature 
range near T... 

A comparison between theory and experiment 
for the longitudinal field variation of the surface 
resistance AR LA) gives the correct amplitude 
of the effect with the same choice @, =2.17, but 
not the correct shape. However, no other value 
of a, (or y) would give a better fit for the data; 
if a, is chosen to fit the shape better, the ampli- 
tude of AR; (Hg) is too small by several orders b 
of magnitude. Thus, a real defect in our present f y 


two-fluid model seems to be indicated. Miller 0 
has shown that nonlocal corrections to the London f y 
equations are important for R, but less impor- a 


tant for X.'° It thus seems likely that a nonlocal f g 
theory for the superconducting electrons is in- 


dicated, at least for AR ; (Ho). be 

The above results seem to indicate that the ta 
excitations above the superconducting state are a 
of several sorts. There are “heavy” electrons ir 
and some “lighter” electrons. Nearly all of the fu 
excited electrons are “heavy” electrons. The So 
“heavy” electrons make the major contribution pa 


to R(0), X(0), specific heat, thermal conductivity, 
and ultrasonic attenuation experiments. The 
“lighter” electrons, being far less numerous, 
have only a small effect on these experiments. 
However, the magnetic field experiments em- 
phasize the contributions of the “light” electrons 
by the weighting factor a,.° Evidence that these 
“light” electrons are connected with the super- 
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conducting state comes from the absence of 

“ight” electrons in the de Haas—van Alphen” 

and cyclotron resonance’? measurements on 

normal tin. At low temperatures (e.g., ¢<0.7) 

one band of “light” electrons is excited (a, =2.17), 
while near T, a second band appears (@,=5). As 
soon as the second band of “lighter” electrons 

is excited, their effect completely dominates 

AX io): The apparent temperature independence 
of a, is compatible with an exponential E vs k for 
the “light” bands of normal electrons. 

It is hoped that a more detailed experimental 
study will be able to decide upon the electronic 
band structure of the superconducting state. 

The authors wish to acknowledge the hospitality 
of the Institute for the Study of Metals, Univer- 
sity of Chicago, where the measurements were 
made. 
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MEASUREMENT OF LOCAL FIELDS AT IMPURITY Fe” ATOMS 
USING THE MOSSBAUER EFFECT 


G. K. Wertheim 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received March 21, 1960; revised manuscript received March 25, 1960) 


The Mossbauer effect! of Fe” has proved to 
bea convenient tool for the investigation of local 
magnetic fields in the transition metals, where 
other methods have not been entirely successful. 
We have investigated the local fields in Fe, Co, 
and Ni, as well as in a number of other sub- 
stances. 

The sources we used were prepared in a num- 
ber of different ways. A source of Co” in me- 
tallic iron was prepared by electroplating from 
acarrier-free CoCl, solution onto a high-purity 
iron foil 0.0025 cm in thickness, and then dif- 
fusing at 950°C in a hydrogen atmosphere. 

Sources in metallic cobalt and nickel were pre- 
pared by evaporating some of the CoCl, solution 
‘0 dryness on a foil of the metal and diffusing as 
before. Each source was then mounted on a loud- 
speaker voice coil driven at constant velocity 

vith a symmetrical sawtooth wave of such mag- 
titude as to produce a 1-mm excursion of the 

‘il. No distinction was made between positive 
ind negative velocity, and counting was continued 
through turn-around. 





As an absorber we have most commonly used a 
0.0025-cm foil of Type 310 stainless steel (25% 
chromium, 20% nickel). This material, although 
paramagnetic, exhibits an unsplit absorption 
line, indicating that the spin correlation time 7 
is sufficiently short to satisfy the relationship 
TA/h<<1, where A is the hyperfine coupling of 
the Fe” nucleus. We have also used diamagnetic 
potassium ferrocyanide which, as expected, also 
exhibits an unsplit absorption line. The salt 
K,Fe(CN),-3H,O was crushed, mixed with a small 
amount of binder, and formed into a uniform 
sheet 0.08 cm thick. This thickness was chosen 
to make the sheet contain the same amount of 
Fe™ per unit area as a metallic iron foil 0.0025 
cm thick. However, the high density of foreign 
matter makes this absorber less attractive than 
the stainless steel. 

Results obtained with the iron source at room 
temperature and the potassium ferrocyanide ab- 
sorber at 78°K, Fig. 1(a), directly exhibit the 
hyperfine structure of Fe”. The three absorption 
lines yield values for the nuclear g factor for the 
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FIG. 1. (a) The hyperfine structure of Fe*’ in me- 
tallic iron obtained with a potassium ferrocyanide ab- 
sorber. (b) The hyperfine structure of Fe’ in me- 
tallic iron obtained with a stainless steel absorber. 
(c) The hyperfine structure of Fe*’ in metallic cobalt 
obtained with a stainless steel absorber. 





first excited state and for the magnetic field at 
the nucleus which are in good agreement with the 


earlier determination by the Argonne group.*® 
The observed effect is reduced by 35% if the ab- 
sorber is at 300°K. Similar results obtained with 
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the stainless steel absorber, for the cases of 
magnetized and unmagnetized iron foil source, 
are shown in Fig. 1(b). The relative intensity o 
the three lines indicates that there is some re- 
manent magnetization in the source which 
strengthens the middle hyperfine line even in the 
absence of an external magnetic field. Figure |( 
shows the results obtained when the iron foil 
source is replaced by a cobalt metal foil. The 
small shift in the location of the hyperfine lines 
indicates that the magnetic field at the Fe™ nu- 
cleus in cobalt metal is 3.1 10° oe, which may 
be compared to the value of 3.3 x 10° oe in iron 
metal. In nickel, at room temperature where th: 
magnetization is 0.95 of the saturation value, the 
corresponding field is 2.6x10° oe. A more de- 
tailed investigation of this question is in progres 

The linewidths in the stainless steel experimer 
are significantly larger than those obtained in th 
original Fe” experiment,’ in which both source 
and absorber were iron metal. The earlier ex- 
periment has been duplicated with our equipment, 
yielding results in agreement with those in the 
literature. Line broadening could arise in the 
stainless steel experiments from a spin correla- 
tion time which is insufficiently short. If this wer 
the dominant effect it should be possible to obtais 
a measure of the spin correlation time in the 
stainless steel through the relation Av>2r7/7A’/’ 
The observed line width, Av, corresponds toa 
correlation time 7 <10-'° sec. On the other hané 
the line broadening might also arise from the 
combining of positive and negative velocities in 
the experimental setup, provided there is a shift 
in the energy of the emitted gamma ray relative 
to resonant energy of the absorber. Such shifts 
may be due to the second order Doppler effects’ 
or due to differences in the environment of the 
iron nuclei in the source and in the absorber. W 
have observed a shift of this nature for Fe™ in 
Ge. 

In addition to our experiments with transition 
metals, we have also observed the Mossbauer ¢& 
fect of Fe” in n-type silicon. The source was 














prepared in a manner similar to that described 
for the cobalt and nickel foils above; as a final 
step, it was etched in hydrofluoric acid to re- 
move any oxide layer which might contain Co” 
a different environment. Consistent results wer 
obtained with both the stainless steel and the 
potassium ferrocyanide absorbers. The data it 
Fig. 2 show that the hyperfine field at an iron 
nucleus occupying a site characteristic of aco 
balt atom in silicon is approximately 3 x10* 0e. 
This value is very much smaller than that found 
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FIG. 2. The hyperfine structure of Fe*’ in n-type 
silicon obtained with a stainless steel absorber. 


for iron in a normal iron site.® The silicon 
source exhibits sufficiently weak hyperfine cou- 
pling so that samples of this type may be useful 
as unsplit sources in the investigation of other 
iron-bearing materials. They have the advantage 


over sources made by incorporating Co™ into 
stainless steel or potassium ferrocyanide that 
self-absorption in the source is negligible due to 
the absence of stable Fe™. 

The author is indebted to W. L. Brown, V. 
Jaccarino, and W. E. Blumberg for discussions 
and suggestions during the course of this work; 
to K. M. Olsen for the high-purity iron, nickel, 
and stainless steel foils used in these investiga- 
tions; and to H. C. Montgomery, G. Feher, and 
J. Eisinger for the loan of numerous pieces of 
equipment. 
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RECOILLESS RAYLEIGH SCATTERING IN SOLIDS 


C. Tzara and R. Barloutaud 
Centre d’ Etudes Nucléaires de Saclay, Gif-sur-Yvette (Seine-et-Oise), France 
(Received March 21, 1960) 


Using the Mossbauer effect, photon sources 
and analyzers extremely selective in energy are 
now available. We study here with such an ana- 
lyzer the recoilless Rayleigh scattering by atoms 
in solids. 

This effect is related to the x-ray diffraction by 
crystais as follows. The interference at the ex- 
act Bragg angles occurs when the scattering is 
elastic with respect to the lattice as a whole, 
that is, without any phonon exchange. Debye and 
Waller have calculated the reduction in intensity 
ofx rays scattered at the Bragg angles ina 
solid at temperature T,! 


3 E ss udu 
R 
epzexr| -— El 4 = } (1) 
2keL4 x%Q ee -1 
where x= 7/6, 6 is the Debye temperature, and 
Ep=(E*/Mc*) (1 - cosé) is the recoil energy 
given to the free atom by a photon of energy E 


scattered at the angle @. 





In the present work, where we detect the elas- 
tic scattering directly by an energy selection 
instead of analyzing a diffraction pattern, the 
factor y7 is the relative number of photons scat- 
tered without energy change. It is clearly the 
same factor which gives the proportion of re- 
coilless y rays in the Méssbauer effect?; in that 
case Ep =E?/2Mc? in Eq. (1). 

In order to measure the factor g7, we have 
studied the Rayleigh scattering for several ma- 
terials: Pt, Al, graphite, and paraffin. The 
23.8-kev photons emitted by Sn’!** were scat- 
tered at 50°+ 5° and absorbed by a Sn’’® foil 
40 mg cm~ thick (almos* completely black for 
the recoilless photons*) (Fig. 1). The scatterers’ 
thicknesses were such that the transmission of 
the y rays was of the order of 10%. 

The Rayleigh-scattered photons are accom- 
panied by inelastically scattered photons (Raman, 
Compton), considerably shifted in energy, so 
that the selective absorption in Sn’!* occurs only 
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FIG. 1. (1) Sn'!** source; (2) scatterer; (3) bismuth 
stopper; (4) 40 mg cm~ Sn!"* foil (71.5 % Sn'!*); (5) 62 
mg cm~? Pd foil absorbing Sn x rays; (6) 1.5 mm 
Nal(T1) scintillator and photomultiplier. 


for a fraction a of all the scattering processes; 
a is extracted from the form factors given by 
Compton and Allison.‘ 

We have measured the relative decrease A of 
counting rate between room temperature 7, = 300°K 
and T,=80°K. The recoilless scattering proportion 
at T, is approximately 

_1 A-e 
" a ne a? 


where f, is the ratio of recoilless emission of the 
Sn'!** source and « is its self-absorption at T,. 
Here f, = 0.32+ 0.015 and «=0.05+0.01. We neg- 
lect the small recoilless emission at 300°K which 
introduces a negligible correction for ?T,° 

The results are given in Table I. The agree- 
ment between the calculated and experimental 
values of gy is reasonably good, especially when 


Table I. Experimental and calculated values of °T, 











1/a A "?, ? calle 
Pt 1.05 0.27 +0.03 0.72 +0.09 0.80 
Al 2.15 0.19 +0.016 0.92+0.09 0.62 
Cc 4.13 0.10 +0.01 0.79 +0.09 0.68 
CH, 5.30 0.020+0.01 0 





we notice that the Debye temperatures are de- 
duced from specific heat measurements rather 
than from x-ray diffraction. 

We have also, using a thin Sn’? foil as a scat- 
terer, observed at low temperature the resonant 
Mossbauer scattering.°® 

This method extends the range of solids which 
can be studied by means of the Méssbauer effect 
or by x-ray diffraction. 

It is a pleasure to acknowledge interesting dis- 
cussions with Dr. Abragam, Dr. Cotton, and 
Dr. Jacrot. 
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STUDY OF THE INTERMEDIATE STATE IN SUPERCONDUCTORS 
USING CERIUM PHOSPHATE GLASS 


Warren DeSorbo 
General Electric Research Laboratory, Schenectady, New York 
(Received March 28, 1960) 


Alers' has recently used the magneto-optic 
Faraday rotation in some cerous nitrate-glycerol 
solutions to observe the intermediate state in a 
superconducting lead alloy. This technique, like 
the bismuth probe? and the superconducting*’* or 
ferromagnetic® powder methods, makes use of 
the diamagnetism of the superconducting state to 
investigate the topography of the superconducting 
and normal domains at the surface of a specimen. 
In all these methods, it is assumed that the local 
magnetic field distribution and variations at the 
surface reflect the amount of normal and super- 
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conducting material in the specimen itself. More 
recently, Alers® has reported large Faraday ro- 
tations in some paramagnetic glasses containing 
cerium meta-phosphate suggesting their use in 
studying the structure of the intermediate state 
in superconductors. 

Utilizing this Faraday effect in similar glass 
prepared by A. Pincus and R. H. Pry of this la- 
boratory, we have obtained some preliminary 
results on resolving the intermediate state in 
various superconducting metals. Included in 
these early studies was the effect of different 
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conditions of metallurgy, strain, purity, crys- 
tallographic orientation, grain size, etc. The 
purpose of this note is to point out that new in- 
formation on the dynamics and time-dependent 
phenomena of the intermediate state not previously 
seen by other techniques may now be observed. 
Because still photographs do not demonstrate the 
full power of this method, motion pictures are 
more useful in recording the kinetic details; e.g., 
wall motions, annealing behavior of frozen-in 
flux, etc. The technique also facilitates studies 
involving both the magnitude and the direction of 
frozen-in moments (not feasible by powder tech- 
niques), Supercooling phenomena and critical 
field determinations. 

The success in the resolution (better than 0.2 
mm) is due in part to the use of very thin glass 
(-0.25 mm thick). Photographs illustrating this 




























FIG. 1. (a) Intermediate-state 
pattern of a tantalum disk (cold- 
worked) at 1.46°K. Happlied 
=460 oersteds. (H-=1060 oer- 
steds.) (b) Intermediate-state 
pattern of a thin tantalum disk 
(cold-rolled) at 1.52°K. Happlied 
=164 oersteds. (H-=1030 oer- 
steds.) 


FIG. 2. (a) Nonequilibrium 
pattern obtained on a rectangular 
plate of single-crystal zone- 
refined tin 10 minutes after trans- 
verse field had been slowly re- 
duced from 222 oersteds to 32 
oersteds at T7=1.46°K. Pattern 
changes with time. (H,=215 
cersteds.) (b) Nonequilibrium 
pattern obtained on a polycrystal- 
line tin disk 21 minutes after 
transverse field had been reduced 
from 325 oersteds to 58 oersteds 
aT=1.40°K. Pattern changes 
with time. (H_- = 264 oersteds. ) 


are presented in Figs. 1 and 2. Light areas are 
normal and the dark regions are superconducting. 
The monochromatic light beam used (5451 A) was 


obtained from a Hg-arc source. The transverse 
magnetic field was provided by a superconducting 
niobium wire-wound solenoid. 

Figure 1(a) shows the intermediate-state struc- 
ture of a polycrystalline tantalum disk (7/8-inch 
diameter 1/16 inch thick), cold-worked by 
compression after machining from a cast ingot. 
Figure 1(b) presents the intermediate-state pat- 
tern of a cold-rolled 7/8-inch diameter disk; the 
disk was cut from a polycrystalline tantalum 
sheet 0.009 inch thick produced by cold-rolling 
of an “as cast” ingot of vacuum arc- melted tan- 
talum. The residual resistance ratios of each of 
these two materials were comparable (i.e., 
R973°K/R4.2°K= 12). These patterns could be 
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altered by subsequent heat treatments which were 
also accompanied by decreases in critical field. 
Irregular patterns similar to Fig. 1(b) have also 
been observed for a cold-rolled sheet of niobium. 
In each of these observations, as the field is 
raised slowly and uniformly from zero to its crit- 
ical value, large and irregular-shaped normal 
domains appear very suddently [see Fig. 1(a) and 
1(b)]. These abrupt jumps of domain walls are 
similar to those observed by Schawlow® on a 
polycrystalline sample of vanadium and on a 
cold-worked sample of rhenium recently reported 
by Schawlow et. al.’ with superconducting powder 
techniques. Experiments to elucidate the nature 
of these boundary motions are now in progress. 
Figure 2(a) shows a domain pattern (nonequi- 
librium) evident in a zone- refined unstrained 
single crystal of tin (rectangular flat plate 1/2 
inch x 1/2 inch x 7/8 inch) about 10 minutes after 
the applied field had been reduced from slightly 
above the critical field (H, = 215 oersteds) down 
to about 32 oersteds at 1.46°K. To reduce the 
field abruptly to zero would result in similar 
patterns, appearing in a shorter time and anneal- 
ing more rapidly; they, therefore, do not photo- 
graph as well. The corrugated pattern in Fig. 2(a) 
is similar, though not identical, to those observed 
by Balashova and Sharvin® on less pure mono- 
crystalline tin spheres. Their observations were 
made at somewhat higher temperatures but were 
reported to have been made under equilibrium 
conditions. Equilibrium patterns of the inter- 
mediate state, in our case, reveal, for applied 
fields greater than about 150 oersteds, a single 
normal domain surrounded by a border area that 
is superconducting. The inner boundary of this 
area seems to move out uniformly with increasing 
applied field. The magnetic lines, presumably, 
pass through the superconducting area in small, 
unobservable channels. The details of these equi- 
librium patterns are not similar to Schawlow 
and Devlin’s® observations on pure tin crystals, 





although the macroscopic anisotropic kinetic ef- 
fects with increase in field observed by us are 
similar to their results. 

Figure 2(b) shows another nonequilibrium pat- 
tern for a polycrystalline (1 mm average grain 
size) disk of tin (7/8-inch diameter x 0.270 inch), 
The sample was made from an ingot supplied by 
the Vulcan Detinning Company and listed as 
“extra pure fine grade” containing 0.0003 % lead 
and 0.0001% iron. These patterns are similar to 
those reported by Schawlow*; however, details 
on their growth and annealing behavior have not 
been observed before. 

Frozen-in flux has also been resolved in zone- 
refined single crystals of lead both in a strained 
and unstrained condition. The patterns are dif- 
ferent in each case. 

Details of the experimental technique and of the 
various results will be presented elsewhere. 

The author has benefitted from discussions 
with many associates, in particular, with R. H. 
Pry for his suggestions and encouragement in the 
use of the glass and with C. P. Bean. He is in- 
debted to D. T. F. Marple and S. Roberts for sug- 
gestions on the optical system. K. T. Aust and 
J. W. Rutter were kind enough to supply him with 
the various single crystals. W. A. Healy capably 
assisted in the experimental work. 
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SPUTTERING THRESHOLDS AND DISPLACEMENT ENERGIES™ 


Robley V. Stuart and Gottfried K. Wehner 
Mechanical Division of General Mills, Minneapolis, Minnesota 
(Received March 28, 1960) 


Sputtering at very low ion energies (below 50 
ev) is of interest in several areas, particularly 
with reference to the question of whether sputter - 
ing threshold energies might be linked to dis- 
placement energies in metals. Other reasons 
for interest are connected with radiation damage 
problems and the sputtering of satellites. The 
major obstacle to work at very low ion energies 
has been the lack of a sensitive and generally 
applicable method of measuring the small num- 
er of sputtered atoms.’ Sensitive electrical 
detectors cannot be employed because sputtered 
atoms are neutral. Surface ionization’ and radio- 
active tracer methods® have been used as sensi- 
tive indicators, but these methods are applicable 
only to a limited number of target materials. 

At very low sputtering yields another difficulty 
encountered is contamination of the target sur - 
face from background impurities in the tube. 

A newly developed spectroscopic method 
eliminates most difficulties and permits sensitive 
relative yield measurements in the low ion en- 
ergy region for any target material. Atoms are 
sputtered from a target in a magnetically con- 
centrated low-pressure mercury or rare gas 
discharge. In mercury the discharge data are: 
gas pressure ~2 microns, discharge current 
between pool-type cathode* and anode ~4 amperes, 
and discharge voltage ~35 volts. The mercury 
pool-type cathode which is very convenient in a 
demountable tube can be retained for noble gas 
studies by freezing the mercury out of the main 
discharge section and supplying a noble gas in- 
stead. The discharge data in argon are: gas 
pressure ~7 microns, discharge current ~4 
amperes, and discharge voltage ~50 volts. The 
sputtered atoms find very favorable excitation 
conditions (spiraling electrons) in the plasma 
and one observes the emission spectrum of the 
larget atoms superimposed on the emission spec- 
trum of the discharge gas. A strong emission 
line of the target material is monitored by use 
ofa monochromator employing a photomultiplier 
detector for intensity measurement. The curve 
of the spectral line intensity plotted as a function 
of bombarding ion energy in most cases closely 
matches the yield curve obtained by absolute 
methods. The sensitivity of this method is such 
that one can readily measure yields done to 10~* 


atom per ion for many metals. Difficulties with 
surface contamination of the target are overcome 
by incorporating in the tube a magnetically oper- 
ated shutter to protect the observation window 
from sputtered deposits and by cleaning the target 
with high ion energy bombardment immediately 
before each measurement. In practice, the shut- 
ter is left closed and the target held at a high 
voltage continuously except for the few seconds 
needed for taking data. 

Yield curves for sputtering by Hg* ions and 
A* ions have been obtained with this method for 
about twenty different metals. Typical curves 
with yields on a log scale are shown in Fig. 1. 
The following features are of interest in these 
results: 

1. Yields at very low ion energy become sen- 
sitive to the anode voltage above a certain criti- 
cal value. This is supposedly caused by the 
formation of doubly charged ions which arrive 
at the target with twice as much kinetic energy 
as singly charged ions. 

2. With the contribution due to doubly charged 
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ions minimized by operating the tube with a low 
anode voltage, one finds that the yields decrease 
with decreasing ion energy more rapidly than 
previously reported. Thus, at least for normal 
ion incidence, sputtering effects at low ion en- 
ergy seem to have been overestimated in the 
past. 

3. By extrapolating the lower part of the yield 
curves, one finds that approximate threshold 
energies for sputtering are in the range 15 to 45 
ev. For argon the product of these thresholds, 
and the energy transfer factor 4m,m,/(m,+mz,)’, 
is with few exceptions (Pd, U, and Th~15 ev) 
close to 25 ev which is the energy generally asso- 
ciated with the displacement threshold for radi- 
ation damage.*° The measurements are presently 


being continued and extended to helium and other 
rare gases. 





‘Work partly supported by the U. S. Air Force Cam- 
bridge Research Center and the Office of Naval Re- 
search. 
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Soviet Phys. JETP 3, 52 (1956)]. 

‘G. K. Wehner, Phys. Rev. 108, 35 (1957). 
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The author’ has shown that a crystal field model 
gives the cubic field splitting parameter a for 
Mn** proportional, to a first approximation, to 
the square of the cubic crystal field. This rea- 
soning may be extended to ions having the half- 
closed configuration 121 +1 and is explained es- 
sentially as follows. Two ions are said to be 
complementary when the patterns of their linear 
crystalline Stark splittings are inverted. This is 
the case for & and /2(2/+1)-*, Now /2/+1 is 
complementary to itself. If its multiplet were 
split, it would show two inverted patterns. But 
since a given multiplet shows only one pattern, 
any multiplet from 121+1 cannot be split bya 
linear action of the crystal field. It can be split 
to the next approximation by the higher order 
perturbation involving squared off-diagonal ele- 
ments of the cubic potential. Thus the result 
may be extended to the /2/+1 configuration. The 
cubic field splitting parameter a would be of the 
form 


a=a(D,) -B(Dg), (1) 


where @ and § are constants, and Dg is a meas- 
ure of the cubic field. In a first approximation 

only the nearest neighbor ions around the S-state 
ion considered may be assumed to contribute to 
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PROPOSAL FOR AN ELECTRON SPIN RESONANCE EXPERIMENT OF S-STATE IONS 
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Dg: Thus one obtains 


vo - (7) (e,e,/R*)(1/336) 






x)).(35 cos*é . -30.cos’6. +3), (2) 
J J J 





where (7*) is the average of r* with respect to 
the radial part of the 7 orbital; R is the inter- 
ionic distance between the S-state ion and neigh- 
bor ions possessing charges e, and e,, respec- 
tively; @; is the azimuth angle of the jth ion with 
respect to the [001] direction, which is chosen 
as the z axis of polar coordinates; and j is 
summed over the nearest ions. It is concluded 
that the parameter a should be proportional to 
the inverse tenth power of R and should be posi- 
tive regardless of the sign of Dg. 

Several electron spin resonance (ESR) experi- 
ments support this result. Baker, Bleaney, and 
Hayes” reported that ESR experiments on Mn™* 
in calcium fluoride yielded 


a =(+0.6+ 0.4) x107* em“. 


On the other hand, using Dg(MgO) - 4.0 xD (CaF, 
and a for MgO (a =+18.6x107* cm), they found 
that from the author’s graph’ of a against Dg thal 
a for CaF, should be ~+1107* em™ which is in 
remarkably good agreement with their value. 
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Matarresse and Kikuchi* did an ESR experi- 
ment on Mn?* in ZnS in which Mn?* is presumed 

to be surrounded by a tetrahedron of four nega- 

tive sulfur ions and the sign of Dg would be, 
therefore, negative. As a result, this system 
provides an important test for the author’s theory. 
They reported a negative sign for a, which seemed 
to be inconclusive since this negative sign was 
determined without measurements of the relative 
intensities of the transitions at low temperature. 
Watkins> made such a measurement at 4°K and 

found that the sign reported was in error. He 
reported that the sign for a should be positive. 

Geschwind® reported that the parameter a for 
Fe’* is positive in both octahedral and tetra- 
hedral O?~ coordination in single crystal of 
yttrium gallium garnet. Moreover, he re- 
examined the spectrum of Fe** in rubidium 
aluminum sulfate where a negative sign for a 
was reported earlier,’ and found that a is posi- 
tive in this case also. These results appear to 
confirm the author’s theory. 

Recently Low and Rosenburger® have observed 
the paramagnetic resonance spectra of Gd°* in 
single crystals of strontium chloride and ob- 
tained 





a(c in their notation) 
= (+39.6+0.1)x107* cm™. 


On the other hand, Low® reported earlier that a 
for Gd** in calcium fluoride is (+185+5) x10 
em“. The ratio of the values of a measured in 
CaF, and SrCl, is about 4.7. The ratio of the 
sr-Cl and Ca-F distances is about 1.28-1.29. 
Assuming that the trivalent gadolinium ion does 
not change this ratio appreciably, they calcu- 
lated the ratio of the values of a in CaF, and 
Cl, obtaining about 11.8. This calculated 
ratio is larger than the experimental one. Thus 


they suspect that there may be a mechanism 
responsible for the splitting being linear in the 
potential. They also point out that another pos- 
sibility is to assume that the Gd-Cl or Gd-F 
distances differ considerably from the Sr-Cl or 
Ca-F distances. 

We wish to call attention to the value of making 
an ESR experiment on S-state ions under high 
hydrostatic pressure which reduces the inter - 
ionic distance R and induces an increase of the 
parameter a. Such an experiment may provide 
an important test for the author’s theory which 
predicts an inverse tenth power dependence of 
aonR. This dependence may be approximately 
replaced by an inverse ten-thirds dependence on 
the volume of a single crystal. If this dependence 
were confirmed, one may calculate the ratio of 
interionic distances between S-state ions, con- 
tained in the host crystal as impurity, and the 
surrounding ions in various crystals. This ratio 
may give important information about the host 
lattice deformation due to impurity ions. On the 
other hand, if an inverse fifth power dependence 
were found it would provide an impoftant clue in 
the investigation of the mechanisms responsible 
for the splitting parameter a. 
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EVIDENCE FOR QUADRUPOLE INTERACTION OF Fe?’™, 
AND INFLUENCE OF CHEMICAL BINDING ON NUCLEAR GAMMA-RAY ENERGY* 


O. C. Kistner and A. W. Sunyar 
Brookhaven National Laboratory, Upton, New York 
(Received March 30, 1960) 


The recoil-free emission and resonant absorp- 
tion’ of the 14.4-kev nuclear gamma ray of Fe*”’, 
has been used to determine the quadrupole coup- 
ling for the 3/2- excited state of Fe™ bound in 
Fe,0,, and to measure an energy shift of this 
nuclear gamma ray which is attributed to effects 
of chemical binding. This effect is corollary to 
the effects of chemical environment on internal 
conversion coefficients? and on electron capture 
disintegration rates.* These measurements also 
yield the value of the internal magnetic field at 
the position of the Fe*’ nucleus when it is bound 
in antiferromagnetic Fe,Q,. 

The M1 emission line of Fe®*” bound in ordinary 
metallic iron is split into six components by the 
magnetic hyperfine interaction. The resonant 
absorption of this emission spectrum by Fe*”’ 
bound in Fe,O, has been examined, as well as 
the much simpler absorption pattern which re- 
sults when the “unsplit” emission line from Fe*”’ 
bound in a stainless steel lattice’ is used. Be- 
cause the interpretation of the two sets of meas- 
urements agree, only the latter measurements 
are presented in this Letter. The former meas- 
urements, in which a Co*”’ source co-plated with 
iron onto 1-mil copper was used, will only 
briefly be remarked upon. 

The ~2-mg/cm? Fe,O, absorber used in these 
measurements contained Fe*”’ enriched to ~30%. 
The source consisted of Co*’ plated onto 0.001- 
inch stainless steel (25% Cr, 20% Ni). After 
plating, this source was annealed for one hour 
at 900°C in a hydrogen atmosphere. All meas- 
urements were made with source and absorber 
at a temperature of 25°C. 

In order to obtain a Doppler shift of the emis- 
sion line, a uniform motion was obtained by 
coupling a pneumatically driven cylinder to 
another cylinder filled with oil, the ports of 
which were connected via a needle valve. A 
wide range of uniform velocities could be selec- 
ted by adjusting this valve. The direction of 
source travel was reversed automatically by 
means of microswitches. Additional micro- 
switches, set to exclude the region of nonuniform 
motion near the travel limits, were used to pro- 
vide gate signals for the counters. The distance 
of travel between the limits of the counting gates 
was 0.973 cm. Source velocity was determined 
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by counting the cycles from a 1000 cps tuning 
fork oscillator during the time between the gate 
limits. 

The 14.4-kev gamma ray was detected witha 
Nal(T1) scintillation counter. The phosphor was 
2 mm thick and 1.5 inches in diameter. The 
counter face was located 5.4 cm above the upper 
limit of vertical travel of the source. The ab- 
sorber was placed 4.7 cm from the counter face, 
A single-channel pulse-height analyzer selected 
the 14.4-kev gamma ray photopeak. The outputs 
of this analyzer and the 1000 cps clock were 
switched between two pairs of scalers so as to 
record counting rate and velocity separately for 
both directions of source motion. 

Figure 1 shows the counting rate (in arbitrary 
units) as a function of source velocity relative to 
the Fe,O, absorber for the stainless steel source. 
Absorption of the “unsplit” emission line at each 
of the six possible absorption energies of Fe™ in 
Fe,O, is evident. The lack of symmetry of the 
absorption pattern about zero velocity shows im- 
mediately that one is not dealing simply with a 
magnetic hyperfine splitting pattern. 

The velocities at which absorption peaks occur 
are given in Table I. Experimental values have 
been corrected by -2.5% to take account of the 
effect of geometry on our velocity scale. These 
absorption line velocities may be fitted precisely 
in terms of an energy level diagram as shown 
in Fig. 2. It is necessary to introduce an energy 
shift AE = AE, + AE, between the center of gravity 
of the absorption lines of Fe*’ in Fe,O, and the 
emission line of Fe®*’ in stainless steel. In addi- 
tion, an energy shift €, of positive or negative 
sign, on the individual m states of the excited 
level is required. This is interpreted as being 
caused by a quadrupole interaction when Fe™ is 
bound in Fe,O,. A least-squares fit to our data 
yields the following splitting parameters (in 
“velocity units”) for the two Fe™ nuclear states 
when Fe™ is bound in Fe,0,: 
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Z,' = 0.611+ 0.005 cm/sec, 
g,' =0.345+ 0.003 cm/sec, 
AE = 0.047 + 0.003 cm/sec, 

€ =0.012+ 0.003 cm/sec. 
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interaction. For a pure quadrupole spectrum, 
in the case of axially symmetric field gradients, 
the shift € for the substates of a spin / state is 
given by® 


2 

€ “Tae pam -KI+1)], 
where q denotes (1/e)(8?V/8z?). All other sym- 
bols have the conventional meanings. When 
Zeeman splitting is combined with quadrupole 
interaction, the precise energy shifts for the 
individual m states due to quadrupole interaction 
depend upon the orientation of the magnetic axis 
relative to the axis of symmetry for the electric 
field gradient, and are in general not the same 
for all m states.” However, for a case such as 
we deal with here, where the quadrupole inter - 
action is small compared to the magnetic inter- 
action, the absolute magnitude of € is the same 
for all m states of the spin 3/2 state within the 
accuracy of our measurements, independent of 
the orientation of the crystal axes. Thus, we 
find € = 4/e7gQ| =0.012+ 0.003 cm/sec (5.75 x10~° 
ev). Since our measurements indicate that the 
shift of magnitude ¢€ is negative for the m =+3/2 
substates and positive for the m =+1/2 substates, 
the product e?¢Q is negative. The absolute value 
is suggestive of a small quadrupole moment for 
Feo?m although this statement cannot be ampli- 
fied further without knowledge of electric field 
gradients in Fe,Q,. 

In our earlier absorption measurements on 
Fe,O, using the emission lines of Fe*’ bound in 
ordinary iron, we find an energy shift upward 
(see Fig. 2) of AE, = 0.04 cm/sec. Within the 
accuracy of our measurements, the splitting 
parameters g,’, g,’, ande remain the same. A 
separate measurement of the absorption by Fe*’ 
bound in stainless steel of the emission lines of 
Fe™ in ordinary iron confirms the implied ex- 
istence of an energy shift downward (see Fig. 2) 
of AE, =0.01 cm/sec in this case. 

The existence of an energy shift AE =(2.26+ 0.15) 
x 10° ev between stainless steel and Fe,O, has 
been definitely established. We will now discuss 
the ways in which differences in chemical environ- 
ment may produce such a shift. (1) Source and 
absorber will in general have different Debye 
temperatures. Since the nucleus in its excited 
state has a slightly greater rest mass than when 
it is in its ground state, the energy difference 
between the nuclear states will decrease from its 
value for an “unbound” system by different 
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amounts because of the difference in the zero- 
point energies. This reduction in energy is lar- 
ger for the substance having the higher Debye 
temperature. (2) In addition, when a lattice is 
at a finite temperature, the energy difference 
between the nuclear states is reduced further by 
virtue of the previously reported temperature 
effect.* Since the magnitude of this reduction 
depends upon the integral of the specific heat of 


the lattice up to the temperature of the substance. 


it will be larger for the substance having the 
lower Debye temperature. We expect this tem- 
perature effect to be small, although it should be 
noted that it acts to reduce the energy shift be- 
tween different substances which results from 
zero-point energy differences. The sign and 
magnitude of the observed shift between stain- 
less steel and Fe,O, would require that Fe,0, 
have a considerably lower Debye temperature 
than stainless steel or ordinary iron. Were this 
the case, the Debye-Waller factor would materi- 
ally depress the recoil-free resonant yield in 
Fe,O,, a result which is not indicated by the data. 
(3) When chemical environment is altered, a 
nuclear isotope shift may result. This effect has 
its origin in the change in the electronic wave 
functions over the region of space occupied by 
the nucleus. s electrons may be expected to con- 
tribute most to this effect. Since s electrons are 
in effect removed in going from Fe in metal to 
Fe in Fe,O,, a smaller charge radius for Fe” 
than for Fe™ in its ground state would produce a 
shift in the observed direction. The direction of 
the observed energy shift AE requires the pre- 
sence of the nuclear isotope shift, since the zero- 
point energy shift is in the opposite direction. 
This nuclear isotope shift is similar to the iso- 

























































meric isotope shift® observed in the optical spec- | 


trum of Hg’. 

We wish to emphasize two additional points 
about recoil-free emission and absorption ex- 
periments when source and absorber are chemi- 
cally different. First, the existence of energy 
shifts introduces asymmetries into the absorp- 
tion pattern which make it essential not to com- 
bine data taken at equal velocities of opposite 
sense. Second, any recoil-free absorption at 
zero velocity is accidental. In general, there- 
fore, a search for this effect by comparing ab- 
sorption at zero velocity with absorption at a 
large relative velocity between source and ab- 
sorber will not yield significant results unless 
source and absorber are identical chemically. I 
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is amusing to note that in the absorption by Fe,O, 

of the emission lines of Fe™ in ordinary iron, 

two such accidental coincidences do indeed occur 

to give substantial absorption at zero velocity. 
We wish to thank many of our colleagues, 

particularly M. Goldhaber and J. Weneser, for 

interesting discussions, and G. K. Wertheim 

for providing us with a sample of the particular 

stainless steel used in his measurements.” 
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STRUCTURE OF NUCLEAR MATTER 


A. W. Overhauser 


Scientific Laboratory, Ford Motor Company, Dearborn, Michigan 


It is usually assumed that the Hartree (or 
Hartree-Fock) ground state of a Fermi gas is the 
familiar sphere of occupied momentum states. 
For a noninteracting gas such is, of course, the 

| ground state; and for an interacting gas, this 
state is, to be sure, an exact solution of the 
Hartree-Fock equations. But to regard it as the 
lowest single-particle state has always been an 
unproved assumption, and indeed a false one. 
Explicit Hartree states of lower energy will be 


; | displayed below for a Fermi gas with attractive 


interactions. This problem is of practical inter- 
est with regard to nuclear matter, and will be 
discussed within that context. It will be shown 
that there are large static density waves in nu- 
clear matter, so that (to an adequate approxima- 
tion) the total nucleon density has the following 
spatial variation: 


P=p,(1+y cosgx)(1+y cosqy)(1+ycosgz). (1) 


In another paper it will be shown that a low- 
energy Hartree-Fock state of a Fermi gas with 
repulsive interactions will have large spin den- 
sity waves. It should not be necessary to em- 
phasize that these new states have not been 

proved to be the ground state in the single-par- 
ticle approximation, although they do have con- 
siderably lower energy than the “normal ground 
state.” However, it seems clear that any attempt 
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to improve upon the single-particle approxima- 
tion, such as the many-body techniques developed 
in recent years, should begin with the lowest 
Hartree-Fock state available. 

It is not our purpose here to present a general 
argument, applicable to all conceivable interac- 
tions, but rather to illustrate the physical struc- 
ture of the new low-energy states within the 
framework of a simple model. We shall treat 
first a one-dimensional problem, and afterwards 
a three-dimensional nucleon gas. Consider a 
large number N of spinless Fermi particles ina 
one-dimensional box of length L. The kinetic 
energy operator will be the usual one. We shall 
assume the total interaction energy of the gas to 
be 

L 
V= -(aL/2n)f p*dx, 


i) 


(2) 


where p(x) is the total particle density and p a 
constant. If the difference between N and N-1 is 
neglected, the Hartree potential U(x) experienced 
by any onc of the particles is 


U(x) = -(BL/N)p(x). (3) 


[The results to be derived are by no means lim- 
ited to the particular interaction (2). For ex- 
ample, any power of p may be used in the inte- 
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grand of (2). However, the requirement that the 
equilibrium density be finite restricts that power 
to values between 1 and 3.] The single-particle 
quantum states will be indexed according to the 
wave vectors allowed by periodic boundary con- 
ditions. The usual momentum states will have 
energy 2wpE pr, where 


wy =3(k/k,)?, (4) 


and Ey is the Fermi energy, h*k,?/2m, with k, 
=nmN/L, the wave number at the Fermi surface. 

Consider now the following single-particle 
wave functions: 


ae rT - 
>, = ze +8, -w,)e “ [erlee,-op)4 . (5a) 


for 0 <k <k,; and 


07 








. ; c] ¥ 


for -k,<k <0. The energy &, appearing para- 
metrically, is 
E, = 2(w, + Wp _¢) . La(w, - pg) +g7|", (6) 

for positive k, and the corresponding expression 
with q=2k, replaced by -q for negative k. The 
foregoing wave functions and energies are the 
solutions of a single particle in a potential, 
U=4gE cos2k,x, (7) 
provided admixture of states with | |>2k, is 
ignored. (These are the wave functions and en- 
ergies which arise in the weak-binding theory of 
energy bands. The energy gap at +k, is 4gF fF.) 
The qualitative feature of the wave functions (5) 
is that their square magnitudes have an oscilla- 
ting component out of phase with the potential 
(7), and give rise to a density wave which can in 
turn yield (7) as the Hartree potential of the 
attractive interactions. States corresponding to 
|k|>k, may be defined similarly, but should not 
be occupied, since their square magnitudes 
oscillate in phase with (7). This is the reason 
why the wave number of the density wave should 
equal the diameter, 2k,, of the Fermi surface. 
(Only those states below the energy gap should 
be filled.) It can easily be proved that this choice 
leads to the lowest total energy. But we must 
first show that a Fermi sea filled with the states 
(5) has lower energy than one filled with plane 
waves. 

The kinetic energy T of a Fermi gas filled with 
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the distorted states (5) can be (tediously) eva- 
luated in a straightforward manner: 


T =5NE ,L1 +3F(g)], (8) 


where 
F(g) =1-(1+g?) +g? ln| g-*+(1+9°?)]. (9) 


The first term of (8) is the kinetic energy of a 
normal gas, whereas the second is the increase 
arising from the assumed distortion of the N- 
particle (simple product) wave function. The 
density of the gas can also be calculated directly, 
and is 


p=(N/L)[1 - G(g) cosqx], (10) 
where 


G( g)= gin[g™*+ (1+ g7?)*]. (11) 


As surmised, the state has a static density wave, 
This result together with (2) determines the in- 
teraction energy: 


V = -4gN(1+4G?). (12) 


If F and G are expanded in powers of g, it can 
be made obvious that V will always decrease 
faster than the increase of T near g=0. Actually, 
the sum of (8) and (12) can be minimized without 
approximation. The optimum g is 


g=1/sinh(4E pr’): (13) 


The resulting change in total energy relative to 
the normal state is 


AW= -NE plcoth(4E p/P) - 1}. (14) 


This result is negative definite. It is also of in- 
terest to note that if (10) and (13) are inserted 
into (3), one obtains a constant term plus the 
oscillating potential (7), the coefficient of the 
latter being in agreement with (13). Consequently, 
the N-particle wave function is a self-consistent 
solution in the Hartree sense (states with |k|>2, 
being ignored). A further observation is that 
within the field of variation employed here, the 
“normal ground state” is an energy maximum. 
The generalization to a three- dimensional gas 
is almost trivial. In order to optimize the effect 
of density waves, the occupied (unperturbed) 
states should fill a cube in momentum space. 
Three density waves will result, according to 
(1), and there will be a large energy gap around 
the surface of the Fermi cube. (The wave num- 
ber q of the waves must equal the edge length of 
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the cube.) The single-particle wave functions 
are merely the products of three factors, each 
having the form (5). To a first approximation 
each of the density waves can be treated inde- 
pendently, and the analysis is identical to that 
derived above. However, employment of an in- 
teraction term such as (2) does not lead to an 
equilibrium density for nuclear matter. 

The energetics of nuclear density waves for a 
model which exhibits nuclear saturation may be 
explored easily with the method of Karplus and 
Watson.! According to their treatment, the bind- 
ing energy per nucleon of the normal state is 
(in Mev) 


B= -(23/1?) + (V,/2n°)[1-(3a/10r?)], (15) 


where 7 is a radius parameter, 7°=p,/p, Pp, being 
the observed nuclear density. The parameters 
V, and a may be estimated from nucleon scatter- 
ing experiments or, for our purposes, may be 
chosen to yield the observed binding energy (~ 16 
Mev) and density (7~1) of nuclear matter. The 
first term of (15) is the kinetic energy per nu- 
cleon, and the second is the potential. The nega- 
tive term in square brackets represents the de- 
crease in attractive interaction associated with 
the mean square relative momentum P of the 
nucleons: 

P~> \k.-k,/?. (16) 
ij tj 
[The physical origin of this term is of course the 
penetration of the repulsive core by nucleons of 
high relative momenta. The potential terms of 
(15) and (2) are equivalent if 8 is allowed to have 
this dependence. Nuclear saturation is thereby 
guaranteed.| The appropriate modification of 
(15) to include the effects of three perpendicular 
density waves is 


3a 
107? 





Bre - (1+3F)+ Pe (ekorrh- age) 


(17) 


where F(g) and G(g) are as defined by (9) and 
(11). For excitations near the Fermi surface, 
the fractional increase in P is one half the frac- 
tional increase in kinetic energy. Equation (17) 
should first be minimized with respect to g. If 
only linear terms in F and G? are employed in 
finding the optimum g, the result is analytic: 


g=1/sinhé, (18) 





where 


6 = (2760n* + 9aV,)/V,(307? - 9a). (19) 


We have computed B and B’ vs n for various sets 
of V, anda. The set, V,~120 Mev, a~1.4, 
yields the observed density and binding energy 
for B’. The extra binding energy (B’ nox - Bmax) 
attributable to the density waves is 4 Mev per 
nucleon. The energy gap, gh?q?/2m, at the 
Fermi surface is 30 Mev. The amplitude y of 
each density wave is 0.6. These quantitative re- 
sults are relatively insensitive to the values of 
V, and a. One may conclude that density waves 
in nuclear matter are of great importance. 

The foregoing arguments have the force of the 
variational theorem. The density waves are so 
large that the variational solution should be im- 
proved appreciably by allowing states with wave 
vector components greater than g to be mixed 
into the single-particle wave functions. Also, 
the optimum value of g will be different from the 
crude estimate (18). Such refinements can only 
enhance the importance of density waves and in- 
crease their contribution to the binding energy. 
In writing (17), we have neglected the 8% in- 
crease in kinetic energy required to fill a cube 
rather than a sphere in momentum space. Most 
of this increase arises from filling the corners 
of the cube. It seems likely that the cube will be 
truncated by four smaller density waves parallel 
to the [111] type directions. The effect on the 
three major density waves should be small, but 
a major fraction of the 8% increase in kinetic 
energy can thereby be eliminated. 

In the limit of weak attractive interactions, the 
Fermi surface will become a many-faced poly- 
hedron, each pair of opposite faces arising from 
a density wave. All single-particle states near 
the Fermi surface will be essentially 100% per- 
turbed compared to the normal state. It follows 
that the low-energy states proposed in this paper 
are not related to the normal state through a 
perturbation expansion in the interaction constant. 
The stability of static density waves arises from 
the fact that the greater binding caused by density 
fluctuations is achieved with less expenditure of 
kinetic energy, when compared with a mere uni- 
form compression of the normal state. This ob- 
servation points out the inadequacy of the Fermi- 
Thomas approximation for such considerations. 

Glassgold, Heckrotte, and Watson? found that 
nuclear matter is unstable with respect to com- 
pressional collective modes. They assumed, of 
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course, that the parentage of the true ground 
state is the normal free particle state. The in- 
terpretation of their result is now clear. Com- 
pressional modes relative to a highly excited 


“ground state” are unstable if they have sufficient 
admixture of low-energy states, similar to those 


considered here. 





'R. Karplus and K. M. Watson, Am. J. Phys. 25, 
641 (1957). 

2a. E. Glassgold, W. Heckrotte, and K. M. Watson, 
Ann. Phys. 6, 1 (1959). 





REGENERATION AND MASS DIFFERENCE OF NEUTRAL K MESONS* 


Francis Muller, Robert W. Birge, William B. Fowler, Robert H. Good, Warner Hirsch, 


Robert P. Matsen, Larry Oswald, Wilson M. Powell, and Howard S. White 
Lawrence Radiation Laboratory, University of California, Berkeley, California 


Oreste Piccioni 
Brookhaven National Laboratory, Upton, New York 
(Received March 29, 1960) 


A very significant feature of the Gell-Mann- 
Pais particle mixture theory’”* is the regenera- 
tion of the K1 from the K2 neutral meson. We 
examine the three possible types of regeneration 


and give the results of an experiment that exhibits 


the expected transformations as demanded by the 
theory. The experiment also allows an estimate 
of the difference between the masses of K1 and 
K2. 

One of the three types of regeneration has been 
described previously*: A plate inserted into a 
parallel beam of K2 particles produces a parallel 
beam of K1 particles. This phenomenon, which 
we will henceforth call transmission-regenera- 
tion, is in striking contrast with other known 
processes whereby a particle transforms into 
another one: a parallel beam of charged pions 
obviously cannot produce a parallel beam of 
neutral pions by interacting with a plate. 


Here we point out another process that typically 


follows from the theory, namely the regeneration 
by diffraction. Because the K° and the K° waves 
are diffracted by a nucleus with different ampli- 
tudes, the diffracted wave contains K1 as well as 
K2 particles. Thus K1 mesons are regenerated 
by a nucleus with a typical diffraction angular 
distribution. 

Regeneration of K1 can also occur by inter- 
action of K2 with single nucleons. The angular 
distribution of this nucleon-regeneration is 
broad, not essentially different from that ob- 
tained in K-nucleon scattering, and therefore 
it is not a crucial consequence of the particle 
mixture theory. 





418 


All three of these components will emerge from 
a plate traversed by a parallel beam of K2’s. 
The angular distribution should permit one to 
recognize each component separately. 

Case* and Good® have shown that the intensity 
of the transmitted component is a very sensitive 
function of the mean life 7, of the K1 and of the 
difference 6m between the masses of K1 and K2. 
The mass difference appears in the final expres- 
sion because of the phase difference it introduces 
between the K1 and the K2 waves, an effect which 
was first noted by Serber® in connection with K° 
production. Moreover, Good pointed out that the 
intensities of both the transmitted and “scattered” 
component (in the forward direction) are pro- 
portional to |f,,°!?, f,,° being the amplitude of 
the regenerated K1, at zero angle, in a K2- 
nucleus collision. Good’s “scattered” component 
must be identified with the diffracted component 
described above. Thus the intensity ratio of the 
transmitted wave to diffracted wave is a function 
only of 5m andt,. We derive here in a more 
concise way the expression for this ratio. 

The computation of the expected transmitted 
and diffracted intensities can be greatly simpli- 
fied by neglecting, from the start, the regenera- 
tion of K2 from K1. As the number of K1’s is 
always less than one thousandth of the number of 
K2’s, this approximation is very good. We con- 
sider then a plane wave of K2 particles, of wave- 
length :, crossing our plate, which contains N 
nuclei per cubic centimeter. If each nucleus 
produces K1’s with a forward amplitude /,,°, a 
infinitesimal thickness dx of the plate at depth * 
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(x in the direction of the incoming K2 beam; 

x-0 and x =L denote the limits of the plate) pro- 
duces a K1 wave amplitude iNAf,,°dx which ar - 
rives at the end of the plate with the amplitude 


, ‘ L-x L 
da, =iNAf .,°dx exp (ie -ik,(L -x) “ion, “&): 
1 


Notice that the K2 wave has traveled to depth 
x before producing the K1 wave, which then 
travels from x to L; u is the collision mean free 
path, which is the same for K1 and K2, because 
both particles are a half-and-half mixture of K° 
and K°; » is the velocity of the particles; y is 
the Lorentz factor; #k, and hk, are the momenta. 
Let us call A =vyt, the decay mean free path of 
the K1’s and introduce the dimensionless quanti- 
ties 1=L/A and 6=(m,-m,)c?/(i/c,). By inte- 
gration with respect to x we obtain, for the trans- 
mitted intensity, 

Alfa? ENE AME | -i6l “1/2 \2 -L/u 
= 12467 e -e e . 


la|? 





(1) 


On the other hand, the nuclei, incoherently from 
each other, regenerate K1’s by diffraction with 

a differential cross section do,, /dw = \|f,,|?. The 
number of diffraction-regenerated K1’s in the 
infinitesimal thickness dx at x, in the forward 
direction, surviving through the thickness L -x is 


dn,\° L-x L 
—' ei 0,2 os 
«(@2) for | Nexp( oT, =) ae 


which is integrated to give 


-L/u 
dw F (2) 


The ratio between (1) and (2) is 
l -1/2 \? 
-e 


o — 
(iz) = If2,°7NA(1 -e , 


R=4NAX" |e * Nia -e (1 +46"). (3) 


To observe these regeneration processes, we 
have inserted a plate in the Berkeley 30-inch 
propane chamber. The chamber was placed in 
a beam of K3 particles which traversed the in- 
strument lengthwise and perpendicular to the 
plate. The experimental setup will be described 
in more detail in a later article. Here we give 
only a brief description. 

A beam of 1.1-Bev/c negative pions impinged 
on a five-foot hydrogen target. The 670-Mev/c 
K° produced in the target travelled a distance 
of 22.5 feet before arriving at the 30-inch pro- 
pane chamber, so that approximately one K2 
crossed the chamber per 10" protons in the 





Bevatron beam. About 200000 pictures were 
taken, half of them with a 1.5-inch iron plate 
to enhance the diffracted wave relative to the 
transmitted wave, and the other half with a 6- 
inch iron plate which yields an intense trans- 
mitted wave. 

We limited the analysis to those two-prong 
events in which the positive prong could be recog- 
nized to be a meson on the basis of ionization 
and momentum. We also required that the decay 
occur within two mean lives from the plate and 
that the primary momentum be equal, within the 
errors, to the beam momentum. The Q(z, 7) 
distribution of these selected events shows a 
marked peak around the expected value of 220 
Mev, which fact proves the regeneration of K1. 
As a further selection, we keep only those events 
for which Q differs by no more than 1.4 standard 
deviations from the peak value. 

By measuring the vector momenta of the two 
prongs, we determine the angle @ between the K1 
and the incident K2 beams, within an error of 
about 2 degrees. 

The angular distributions are shown in Fig. 1. 
The diffraction curve has been computed with an 
optical model method’ using the known cross 
sections for K* and K~ with protons and nuclei.® 
The curve is quite close to the black-sphere dis- 
tribution. Figure 1 clearly shows a diffraction 
component and a superimposed transmission 
peak. Most of the transmission peak is confined 
to angles smaller than 2.5 degrees (cos >0.999), 
which is just what we expect from an infinitely 
narrow peak measured with our errors. The 
mere presence of such a large transmission peak 
is a proof that the mass difference is smaller 
than, say, 5f/r,. 

Referring to the data for the 6-inch plate, 29 
events occur in the interval for cos@ between 
0.998 (3.5°) and 1, which should contain the total 
number T of the transmitted K1’s and a part D 
of the diffracted ones. Knowing the diffraction 
angular distribution we can compute D from the 
31 particles in the 0.980 to 0.998 interval, as- 
suming that that region contains only diffracted 
K1’s. A nucleonic background will actually be 
present in this interval, which we ignore for the 
moment. We thus obtain D=12, which gives 
T=17. For comparison with formula (3), we note 
that (dn, /dw) =1.18(D/w), where the factor 1.18 
is equal to the ratio of the intensity at zero de- 
grees to the average intensity within w, w being 
the solid angle in the peak interval, that is, 
2n X0.002. It is convenient to compare 7T/1.18D 
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FIG. 1. Histograms of number of K1 decay events 


per 0.001 interval of cos@ (6 is the angle between the 
direction of the primary K2 beam and the regenerated 
Ki). (a) Data for the 1.5-inch plate; (b) data for the 
6-inch plate; (c) combined data for the two plates. 
The curves are diffraction angular distributions nor- 
malized in the 0.980 to 0.998 interval for cos @. 


with R/w, which is plotted in Fig. 2 versus the 
mass difference, in units of i/7,. We obtain 
T/1.18D =1.2+0.53, which gives 5m =0.8579°4,. 
With a probability of 95%, 6m<1.4. 

In the same way, we find from the thin plate 
alone 6m =0 (6m<4.5 with 95% probability) and 
from the combination of thin and thick plate 
5m =0.8579-3, (6m <1.5 with 95% probability). 

If we correct for nucleonic or any other back- 
ground, we would obtain a larger value for R, 
hence a smaller value for 5m. For instance, 
assuming a uniform background from cosé@ = 0.96 
to cos@ =1, the compounded data for both plates 
yield 5m =0 (6m<1.1 with 95% probability). In 
view of a remark by Okun’ and Pontecorvo,° this 
result indicates that decay rates for AS =2 are 
10° times slower than for AS =1. 
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FIG. 2. Calculated intensity of the forward regen- 
erated K1’s versus the K1-K2 mass difference. The 
ordinate R/w is the ratio of the transmission-regen- 
erated K1’s to the diffraction-regenerated K1’s in the 
interval cos 6>0.998. 
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COSMIC RADIATION INTENSITY DECREASES OBSERVED AT THE EARTH 
AND IN THE NEARBY PLANETARY MEDIUM* 


C. Y. Fan, P. Meyer, and J. A. Simpson 
Enrico Fermi Institute for Nuclear Studies, University of Chicago, Chicago, Illinois 
(Received March 29, 1960) 





This note reports preliminary experiments presence of the earth and its geomagnetic field. 
which prove that existing theories for the modu- On the other hand are models which depend upon 
lation mechanism responsible for rapid decreases the solid earth and its permanent magnetic field 
of primary cosmic-ray intensity cannot invoke for creating the decrease of cosmic-ray intensity; 
the presence of the earth or its magnetic field. namely, geocentric models.*® To decide between 

Forbush first noted sudden intensity decreases these two classes of hypotheses, we ask the 
in ionization chambers located deep within the question: How far into the interplanetary medium 
atmosphere which followed in the order of 20 to is the full decrease of galactic cosmic-ray 
40 hours some large solar flares.‘ These in- intensity observed during a Forbush decrease ? 
tensity decreases arose from changes in the Is the pre-existent radiation intensity found at 
secondary radiation produced by primary cosmic distances beyond which the geomagnetic field 
rays and were frequently accompanied by sub- could be invoked to account for the Forbush 
stantial geomagnetic field disturbances. It was decrease, or is the intensity also reduced in the 
later shown experimentally that this kind of nearby interplanetary medium as it is at the 
world-wide decrease was a property of the pri- earth? For all existing models which invoke the 
mary cosmic radiation observed at the earth, presence of the geomagnetic field, a suitable 
and arose neither as a consequence of the geo- limit for the effective extent of the modulating 
magnetic storm nor by phenomena which might region is 6 to 10 Rg, where Rez is the radius of 
change temporarily the magnetic cutoff rigidities, the earth. 
such as ring currents around the earth.? The To answer these questions, we have performed 
changes in the primary spectrum during times a direct experiment using a cosmic-ray detector 
of drepressed intensity~namely, a larger de- carried by the Explorer VI satellite launched 
crease for low magnetic rigidity particles than August 7, 1959 in an orbit which extended beyond 
for high rigidity particles— suggested that the the region of trapped Van Allen radiation and 
modulation was due to interplanetary magnetic reached a range in excess of 7.5R,. In addition 
fields in the vicinity of the earth.? For a typical to the satellite observations, we measured 
event, the high rate for reduction of prevailing simultaneously the changes in the nucleonic 
cosmic radiation intensity at the earth also placed component at the earth over a wide range of geo- 
bounds on the requirements for magnetic field magnetic cutoff rigidities. For this purpose we 
intensity and scale size of suitable modulating used neutron intensity monitors extending from 
electromagnetic fields, i.e., intense fields of the geomagnetic equator to high latitudes. 
large scale size, or vice versa. The cosmic-ray detector in the satellite was 

Several hypotheses have been advanced for a triple-coincidence, proportional counter sys- 
this solar-produced modulation of galactic cosmic tem which measured protons in excess of 75-Mev 
tadiation. On the one hand are models which energy, or electrons in excess of 13-Mev energy. 
may be described roughly as heliocentric; namely, The triple coincidence detector does not respond 
vhere the cosmic radiation intensity in a sub- to bremsstrahlung, such as from electrons trap- 
stantial volume of the inner solar system is ped in the geomagnetic field. The satellite orbit 
reduced by either disordered, * or ordered mag- was an ellipse of apogee approximately 48 800 
etic fields*» 5 of solar origin independent of the km and perigee approximately 6600 km, with 
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its orbital plane inclined by about 46° to the geo- 


graphic equator. A description of the instrument, 
its performance, and the satellite and its motions 


will appear elsewhere.” Most of the data rec- 
orded at the earth during the two months when 
the telemetry was operative covered the first 
four weeks after the launching. During this 
time, one outstanding sharp decrease of cosmic- 
ray intensity and two smaller decreases were 
observed in the neutron intensity monitors at 
the earth. The results we report here are for 
this period. Only satellite data in the range 
35000 to 48 800 km were used for analysis. The 
continuous record of the neutron intensity mon- 
itor at Climax, Colorado (located at a geo- 
magnetic rigidity cutoff of ~2.5 Bv) was chosen 
from among the network detectors to represent 
the changes at the earth. Only neutron intensity 
data for time intervals which overlapped the 
satellite triple-coincidence data were used as 
shown in Fig. 1. The statistical errors of the 
satellite data are less than 1% for all points 
used in Fig. 1. All neutron data have smaller 
errors. From Fig. 1 the over-all correlation 
between the detectors at 35000 to 48000 km 

and at the earth’s surface is obvious. Only the 
large decrease on August 20 provides data of 
sufficient accuracy to compute the ratio of the 


___ AVERAGE TRIPLE COINCIDENCE RATE IN 
SATELLITE FOR DISTANCES >35x1I0> KM 


__" ___ COUNTING RATE OF NEUTRON MONITOR 
510}— AT CLIMAX lero 
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FIG. 1. Cosmic-radiation intensity as a function 
of time showing on August 20, 1959 a Forbush-type 
intensity decrease of ~15% within five hours. The 
nucleonic component monitor provides a measure of 
primary intensity changes at the earth. The light- 
dashed and solid lines are used to connect the periods. 
for which satellite data are available. 
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FIG. 2. The relative changes in cosmic-ray inten- 
sity at the satellite and at the earth for the sharp in- 
tensity decrease of August 20 shown in Fig. 1. The 
correlation coefficient is 0.96. 


relative change of intensity at the satellite and 
at the earth as shown in Fig. 2. The ratio is: 


Relative change in satellite detector 
: : 21%. 
Rel. change in nucleonic component detector 





The correlation coefficient for the two sets of 
data in Fig. 2 is 0.96. We now compare the 
amplitude of the Forbush decrease at great dis- 
tances from the earth with the amplitude in the 
primary radiation at the top of the atmosphere. 
Therefore, we have to relate the intensity changes 
of secondary neutrons as observed at Climax to 
the changes of primary intensity outside the 
atmosphere above Climax. For this purpose we 
may use neutron intensity measurements at 
balloon altitudes which have been carried out 
for typical Forbush decreases over a number of 
years by our laboratory. For the present phase 
of the solar cycle, and cosmic-ray spectrum, 
this factor is approximately 1.8. Therefore, 
within an experimental error of approximately 
10%, the magnitude of the sharp decrease of 
primary cosmic radiation was the same at the 
earth as it was out to distances as great as 
7.5R,. Thus the decrease of August 20, 1959 
means that more than 15% of all the cosmic-ray 
particles were removed within 5 hours by the 
modulation mechanism. Using neutron intensity 
monitor data we have shown that this event is 
worldwide, and that the spectrum has undergone 
changes which are typical of other Forbush 
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decreases Studied at this phase of the solar cycle. 

We believe this experiment excludes existing 
hypotheses which invoke the presence of the 
earth and the geomagnetic field to account for 
the Forbush decreases of cosmic-ray intensity. 
Thus, geocentric models such as solar plasmas 
containing disordered magnetic fields overlying 
the geomagnetic field, or the roughing-up of the 
geomagnetic field by a solar wind do not account 
for our observations. On the other hand, we 
know that these sharp intensity decreases at the 
earth are the delayed consequences of large solar 
flare outbursts on the visible side of the sun. 
Forbush decreases are observed at the earth 
only for restricted angular positions of the solar 
events relative to the position of the earth. 
Therefore, the modulation region may be con- 
fined to a volume of space surrounding the 
Sun-Earth line. This restricted heliocentric 
system of fields is different, at least in scale 
size, from the mechanism required for the mod- 
ulation of cosmic radiation intensity over the 
ll-year solar activity cycle.* We do not here 
go into further details of suitable models. 

The space probe Pioneer V (launched March 
11, 1960) carrying equipment identical with the 
Explorer VI satellite may soon provide further 
tests of the scale size for the Forbush intensity 
decrease at distances >> 8Rp. 

Figure 1 also shows that the radiation meas- 
ured by the satellite undergoes a stepwise in- 
tensity decrease following each of the two de- 
creases which is not apparent in the relativistic 
energy range of the spectrum as represented by 
the neutron intensity monitors on the earth. 
This may indicate the persistent exclusion of 
some low-energy particles by the mechanism 
producing the Forbush decreases and is being 
investigated. 
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ANALYSIS OF THE EXPERIMENTAL 7’+ DECAY SPECTRUM AS A TEST OF THE AT=3 RULE* 


S. Bjorklund, E. L. Koller, and S. Taylor 
Department of Physics, Stevens Institute of Technology, Hoboken, New Jersey 
(Received March 23, 1960) 


A further test of the isotopic spin selection 
rule AT =$ has been recently proposed by Wein- 
berg.’ For the form of the energy spectrum of 
the odd pion in 7 or tT’ decay, Weinberg writes: 


T 3\_ 3m. / 1 4 aT. a 

p Me Q \L+aQ/3m,, "me, * max my P 

(1) 
where T, is the kinetic energy of the odd pion, Q 
is the total kinetic energy of the three pions, 
my is the mass of the K meson, and x,,,, is the 
maximum value of the x coordinate of the Dalitz 
plot. Weinberg shows that if the decay follows 
the AT=4 rule, then the value of a for the 7’ de- 
cay is uniquely determined by the value of a for 
the t decay. From a recent study of 959 7 events? 
Weinberg obtains a,~=4.9 and predicts a_, = -2a, 
= 9.8. 

However, this value for a_, is not consistent 
with taking only a linear term in T, in Eq. (1), as 
p would become negative for large T,. Hence, 
the smallest value of a, consistent with p re- 
maining everywhere positive, a-,=-9.3, has been 
used in this communication as the predicted 
value. This value is well within the errors, 
since the intrinsic uncertainty in Weinberg’s pre- 
diction, due to neglect of quadratic terms in 
Eq. (1), is stated to be about 20%, while the ex- 
perimental error is also about 20%. 

Seventy-two 7’ decays, previously reported by 
the Columbia Emulsions Group,® have been re- 
analyzed to measure the value of a;,. The likeli- 
hood function obtained from Eq. (1), after nor- 
malization, is 








72 
P= [J p; 
i=1 
K 2 72 a oP 
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rT’ i=1 K K 


where K is a constant which does not depend on 
a,,. The relativistic boundary of the Dalitz plot 
was used in the normalization. A plot of the 
likelihood function is shown in Fig. 1. The maxi- 
mum likelihood value for a_, is -7.1, as com- 
pared to the value of -9.3 predicted by Weinberg’s 
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bability that our seventy-two events turn out the way 
they did, has been plotted as a function of the param- 
eter ay 


analysis. At a,,=-9.3, the likelihood function has 
dropped by a factor of about 6. At a value of 
a,-=4.6 (minus one half the predicted value of 
-9.3 and hence corresponding to the 7 value), the 
likelihood function is down by a factor of about 
270. 

Six-division x? tests have also been carried out 
for various values of a,,, yielding the results 
given in Table I. The sensitivity of the data to 
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Table 1. Results of the x? tests of the data for var- 
jous assumed values of the parameter a,:. 





2 





a, X Probability Remarks 
4.6 13.6 0.017 T value 
0.0 © ©68. 8® 0.12 
-1.1 3.8 0.57 maximum likelihood 
value 
-9.3 8.2 0.14 predicted value 





“This value of x? differs from that in reference 3, 
since here the relativistic boundary of the Dalitz plot 
was used, while reference 3 used the nonrelativistic 
approximation. 


errors in measuring the kinetic energies has 

been tested by recalculating the maximum likeli- 
hood value of a, using the upper and lower error 
limits of the most sensitive event of the data, 
which also happens to have large errors. This 
event has T,=49.0*4:3 Mev. The maximum likeli- 
hood values so obtained are a,-= -6.8 and a, 


= -7.4 for the upper and lower limits, respect- 
ively. 

The measured value a,,=-7.1 is in satisfactory 
agreement with the predicted value of -9.3. 
Hence, the present data, as tested by Weinberg’s | 
analysis, are consistent with a AT=4 rule. It is 
unlikely that a,, is as large as a;, although the 
possibility cannot be completely excluded by the 
present data. 

We thank Mr. T. Faith for his aid in carrying 
out the calculations. 
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AUGER EFFECT IN MESONIC ATOMS 


Richard A. Ferrell* 
University of Maryland, College Park, Maryland and CERN, Geneva, Switzerland 
(Received March 21, 1960) 


The Auger effect is important in numerous ex- 
periments dealing with the capture of negative 
particles by nuclei, because of the dependence of 
the interpretation of such experiments on the de- 
tails of the atomic cascade process. It is often 
essential to know, in the case of a strongly in- 
teracting negative particle, in which state it un- 
dergoes interaction with the nucleus. The answer 
to this question depends in turn upon the transi- 
tion rates out of the various states. Therefore 
itis necessary to have a reliable theory of the 
Auger transition rate. Unfortunately, although 
the usual perturbation theory of the Auger effect 
is straightforward, it has been put into doubt by 
the serious discrepancies found by Stearns and 
Stearns‘ for the x-ray yields in .- and 7-mesonic 
atoms. The purpose of the present note is to 
point out an alternative way of computing the 
Auger effect which makes it clear that the usual 
perturbation treatment is absolutely reliable. In 
fact, we demonstrate a quantitative relationship 
between the Auger rate on the one hand and the 





photoelectric cross section on the other. From 
this relationship, it is possible to calculate the 
Auger effect directly from the empirical x-ray 
absorption data. As the Auger rate is too small 
and no other significant radiationless process 
exists, it follows that the effect found by Stearns 
and Stearns must have its origin outside the tar- 
get. This conclusion is supported by noting that 
the discrepancies for the u-K, w-L, and 1-L 
series, although usually considered separately, 
are really the same and are a function of only 
the photon energy, and independent of the target.’ 
For example, the targets of boron, fluorine, and 
sodium give, respectively, u-K, 7-L, and p-L 
X rays all with photon energies within 5 kev of 

47 kev, and all with a yield of only one-half of 
the expected yield. 

The Auger rates for the u-K and yp -L series 
have been computed by Burbidge and de Borde.* 
These rates are too small by factors of 300 and 
of 30, respectively, to explain the drop in x-ray 
yield found by Stearns and Stearns. It might be 


425 




















































VOLUME 4, NUMBER 8 


PHYSICAL REVIEW LETTERS 





APRIL 15, 19% 









supposed that somehow the Auger rate calcula- 
tion could be in error to this extent (although 
this would be surprising because of its straight- 
forwardness).* For this reason we wish to em- 
phasize here the reliability of the theoretical 
Auger rate, by pointing out the following entirely 
new way of making the calculation. 

The combination of the central Coulomb field 
of strength Z-1 centered at the nucleus and the 
perturbing dipole field of the meson is equivalent 
to a new central Coulomb field shifted in position 
to the charge center of the nucleon-meson sys- 
tem. Since the charge center oscillates with the 
transition frequency w, it is convenient to work 
in the accelerated frame of reference attached 
to it. In this frame an oscillating inertia, or 
gravitational, force acts on the atomic electrons. 
The effect of this force is the same as that of a 
uniform fluctuating electric field, allowance 
being made for the charge-to-mass ratio. Thus 
the Auger rate is equal to the rate of photoelec- 
tric ionization, when the atom is subjected to an 
“equivalent photon flux” given by 


c 


l==—>— mY (2-1)? DIF r 
2thw \ e f fi 


=(Z-1)% 07," w,, (1) 
where we made use of the standard formula for 
the radiation rate wp, and have introduced the 
Thompson cross section o7 = (87/3) (e?/mc?)?. F 
is the meson coordinate and the other symbols 
have their conventional significance. Since the 
ionization rate is given by multiplying the photon 
flux by the photoelectric cross section 07_;, we 
find for the ratio of the Auger to radiative rates 


/ O76 (2) 


- _-1)72 
a a tale = te 


It should be emphasized that despite this equi- 
valence of the Auger rate to the photoeffect, the 
Auger effect is not in fact due to the conversion 
of photons, and does not take place in the outer 
radiation zone of the mesonic atom. It instead 
takes place in the near zone and is due only to 
the quasi-static electric field. Indeed, Eq. (2) is 
valid only when retardation is negligible; but 
this is true in the cases of interest. Of course, 
Eq. (2) is also based on treating the meson as a 
dipole. Since the meson orbits are one to two 
orders of magnitude smaller than that of the 
electron, this is a justifiable approximation. 

One can use Eq. (2) to calculate the Auger rate 
directly from x-ray absorption data. Consider 
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the L x ray emitted by a » meson captured in 
fluorine, of energy 32 kev. The Z-1 atom, oxy- 
gen, has an interpolated mass absorption coeffi- 
cient® of 0.302 at this photon energy, correspon. 
ing to a photon cross section of 8.06 <10-™ cm’, 
Subtracting from this eight times the Compton 
cross section® of 0.601 x10~** cm? we find 


=3.25x107"4 cm’, 


°Z-1 
or from Eq. (2) 
w ,/w p= 0.075. (3) 


This can be compared with the value of Burbidge 
and de Borde’ of 0.083. But the latter authors 
did not allow for the screening of one K electron 
by the other. The screening is approximately 
taken into account® by further reducing the effec- 
tive nuclear charge by 5/16, so that the electron 
density at the nucleus must be reduced by the 
factor 


[(Z-1 - 5/16)/(Z-1)]8 = 1-15/16(Z -1) = 0.88. 


We thus obtain the theoretical value of 0.073, in 
good agreement with Eq. (3). As this example 
demonstrates, Eq. (2), by making use of the 
empirical x-ray absorption data to avoid the te- 
dious calculations with the free-electron Coulom 
wave functions, may have some practical advan- 
tage. Equation (2) has the further advantage of 
automatically including screening and the contri- 
bution of the higher electron shells, when these 
are important. 

To continue with the example of fluorine, the 
u-L yield found by Stearns and Stearns was only 
(24+2)%. The u-K series yield for the same 
atom indicates that approximately 20% of the 
muons jump directly to the ground state from the 
n=3 and higher shells, so that the L-yield would 
be expected to be about 80%.° To reduce this 
figure to 24%, the Auger rate would have to be 
31 times greater than given by Eq. (3). In view 
of Eq. (2), this would require a photoelectric 
cross section greater by the same factor, which 
is, of course, excluded. On the other hand, as 
emphasized by Bernstein and Wu,’° processes 
involving electrons of neighboring atoms make 
no significant contribution, since the only essen- 
tial feature of the initial electron distribution 
which enters is the density at the nucleus. The 
further possibility of low-energy radiationless 
transitions between levels belonging to the same 
quantum number n=3, as suggested by Krall and 
Gerjuoy™ for n=2, can be ruled out as there is 
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n0 metastable state available.” In any case, 
from the existence of the K radiation the mesons 
are known to have passed through the n =2 level. 
Hence, one can conclude that the L radiation must 
have been emitted but somehow eluded detection. 
To test this hypothesis, it is natural to assume 
that the discrepancies found with the u-K, yu-L, 
and a-L yields, although generally discussed 
separately, are in fact the same thing. In order 
to compare the various yield curves, it is neces- 
sary to plot them vs energy rather than vs atomic 
number.’* Figure 1 shows the data of Stearns and 
Stearns: and of Stearns, Stearns, and Leipuner, “* 
plotted in this manner. As a rough approximation 
all three sets of data, taken with the same appa- 
ratus, define the same curve. This indicates a 
yield independent of the target and depending 
only upon the quantum energy, or, “by defini- 
tion,” an instrumental effect. A reasonable fit is 
given by a linear fall-off of detector efficiency 
with decreasing energy, beginning at about 90 
kev. Taking such a linear efficiency and dividing 
itinto the 7-K yields found by Stearns and 
Stearns*® changes these into a decreasing mono- 
tonic function of atomic number, as expected if 
only nuclear absorption competes with radiation. 
Furthermore, dividing this linear efficiency 
function into the roughly constant background 
found below 90 kev converts it into a brems- 
strahlung-type spectrum similar to that found for 
the background above this energy.*® 
It is well known that the decreased efficiency 
cannot be attributed to the sodium iodide crys- 
tals used for detecting the x rays. Neither can 
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FIG. 1. u-K, u-L, and 1-L x-ray yields vs photon 
energy, iw. The three sets of data, taken with the 
same experimental apparatus but with different target 
materials, can be accounted for by the single linear 
efficiency curve shown. 





such a linear efficiency curve be attributed to 
some absorber overlooked between the targets 
and the detector. Any x-ray absorber which 
gives a break at 90 kev gives too strong an ab- 
sorption at the lower energies. But it is some- 
times possible for an inefficiency to be intro- 
duced by the electronic circuitry for small pulse 
heights.*” Finally, note should be taken of the 
decrease in yield reported in reference 14 for 
the higher cyclotron intensities. It would be de- 
sirable to have further experimental studies of 
this effect, especially in the crucial energy range 
of less than 90 kev. 
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in the experiments to guard against such a possibility 
[M. B. Stearns (private communication)]. It should 
also be noted that the Rochester group independently 
found a similar decrease in1-L x-ray yield, but 





with a somewhat weaker fall-off at the lower energies 
{M. Camac, M. L. Halbert, and J. B. Platt, Phys, 
Rev. 99, 905 (1955)]. One can only conclude that 
similar experimental conditions must have prevailed, 





UPPER LIMIT FOR PRODUCTION OF =" HYPERFRAGMENTS BY 
K~ CAPTURE IN DEUTERIUM* 


Orin Dahl, Nahmin Horwitz, Donald Miller, and Joseph Murray 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received March 21, 1960) 


Recently, Gondolfi et al. reported evidence for 


the existence of a bound state of the =" system.’ 
Pais and Treiman had already suggested that such 
a system might be produced when K” mesons were 


captured in deuterium.” They assumed that the 

=n would form a deuteron-like structure bound 
in either the ‘S or °S state and used an impulse- 
type model to estimate the relative rates for the 
processes 


K~+d-Z"n+n", (1) 


reactions (1) and (2), we are examining 2100 
=~(n) productions obtained during two exposures 
of the Alvarez 15-in. deuterium chamber to the 
separated 450-Mev/c K~ beam. Of the > (n) 
productions, 85+ 5% result from K absorptions 
at rest. Thus far we have analyzed in detail 227 
events in which the = (m) came to rest and was 
captured via the reactions* 


D"(n)+d-A (or 5°)+n+n(+n). 
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(2) For the three cases giving best fits when these 
events are interpreted as examples of reaction 
(1), the probabilities for exceeding the observed 
values of xy? were 30%, 10%, and1%. From this 
result it is concluded that the fraction of © pro- 
ductions leading to formation of the bound state 
is <1%. 

Day, Snow, and Sucher have argued that it is 
extremely likely that K” mesons, when stopped 
in liquid H, or D,, are captured from high-lying 
S orbitals.® If their prediction is correct and if 
the existence of the >°m is confirmed, the present 


K°> +d=2"+n+n", 


for both odd and even K-> parities. Day and 
Snow have pointed out that the strong S-wave K - 
nucleon interaction will dominate the capture pro- 
cess from either S or P atomic orbitals and have 
extended the calculations of Pais and Treiman to 
the latter case.* Typical values calculated for 
the fraction of all = productions leading to the 
formation of the bound state are summarized in 
Table I. 

In order to determine the relative rates for 


Table I. Probability for formation of a £~ with binding energy Ef} when a K~ meson is captured in deuterium 
via the S-wave K-nucleon interaction. 








=n 
K- bound S-orbit capture* P-orbit capture? 
parity state Eb=0.1 Mev Ep=0.5 Mev Ep=0.1 Mev Ep=0.5 Mev 























0.20 0.37 0.03 0.06 


0 





forbidden 0 











8See reference 2. 
bsee reference 3. 
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result can be consistent with the calculations of 
Pais and Treiman only if (a) the binding energy 
of the =m is exceedingly small (£5 kev) or (b) 
the K - parity is odd so that reaction (1) is for- 
pidden by conservation of angular momentum and 
parity. ’ 

In the absence of any special mechanism, the 
usual de-excitation processes for the K -d atom 
(collisional Auger effect and radiation) will lead 
to population of the 2P atomic orbital from which 
a large fraction of K ~’s will undergo direct nu- 
clear capture. Since the hyperfragment rates 
are suppressed by a factor ~10 (see Table I), no 
strong conclusion could then be drawn from the 
present experiment. 

The following considerations lead us to believe 
that the = "(z) capture events are more favorable 
than decay events for correct identification of 
hyperfragments: 

(a) Ross has examined the interactions of D~ 
hyperons produced by absorption of stopped K~ 
mesons in hydrogen.® He finds that more than 
95% of the > captures occur at rest. Therefore, 
it may be assumed that the present analysis is 
effectively limited to events in which the = (n) 
does not interact (or decay) before traveling its 
full range. If reaction (1) is produced by a 
stopped K~, the =~» will be collinear with the 7* 
and have a range of 2.43+0.05 mm. Since the 
diameter of such a heavily ionizing stub is ~ 0.3 
mm, the usefulness of the collinearity require- 
ment would decrease rapidly if the shorter path 
lengths of the decay events were included in the 
analysis. 

(b) The distribution in the =" -1* included angle 
for reaction (2) peaks at 180°, and many events 
inwhich the > decays within 2.5 mm of the pro- 
duction vertex would be indistinguishable from 
those of (1). Therefore, the present requirement 
of unique range provides a strong additional cri- 
terion for correct identification. 

(c) No corrections are necessary for (i) events 
in which the 5" (n) range is poorly determined 
because of emission of a decay pion at a small 
angle or (ii) collinear £~’s produced by aK” ab- 
sorption in the hydrogen contamination of the 
bubble chamber. 

The film coordinates of all tracks associated 
with the 5“ -capture events were punched onto 
IBM cards with the LRL measuring projector and 
the tracks were spatially reconstructed by means 
of suitable IBM programs. The distribution in 
Tange and included angle for the reconstructed 
events (with included angle greater than 140 deg 
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FIG. 1. Measured range vs £~() -7* included angle 
for reactions (1) and (2). The angle and range re- 
quired for a 2” hyperfragment are designated by P. 
Typical measurement errors are indicated for the six 
closest events. 


Ig0 


and => range between 1 and 6 mm) is shown in 
Fig. 1. There is no evidence for increased con- 
centration of events in the region expected for 
hyperfragments. As a check on the machine 
measurements, all events were re-examined on 
the scanning projector and obviously noncollinear 
events rejected. For the remainder, the © (n) 
was assumed to be collinear with the 7* so that 
the measured values of the 7* dip and azimuthal 
angles (determined to + 0.5 deg) could be used to 
calculate the lengths expected on the film for a 
range of 2.43 mm. Comparison of calculated 
lengths with direct measurements on the film 
gave results consistent with the x? estimates de- 
duced from the machine measurements. 

The significance of the present measurement 
depends upon the probability P, that a Dn will 
come to rest and be captured. This quantity may 
be written 


D 


P=exp[-0 (7, +7 a )], (3) 
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where I, is the decay rate for the bound =, 7, 
is the proper time from production to the time 
the ="n has slowed down to a velocity v/c = 0.05 
(the approximate lower limit for the validity of 
the usual range-momentum relation), 7, is the 
time from v/c = 0.05 to initial capture by a D, 
molecule, and T,P is the time from molecular 
capture to nuclear absorption. 

Snow has shown that if charge independence is 
valid, the = n will probably be bound by less 
than 0.5 Mev.” Dalitz and Liu have argued that 
for such loosely bound systems, the major modi- 
fications in decay rates should arise from the 
change in phase space available and from the 
operation of the Pauli principle in the final state.® 
Using Eq. (3.3) given in their paper, we have 
estimated the magnitude of these effects. Since 
the S- and P-channel decay amplitudes for the 
free = are not known, results for [ are given 
in Table II for the extreme cases. 

Integration of the range-momentum relation 
yields T,=0.99x10"*° sec. Wightman has esti- 
mated T, for particles of various masses; ° a 
slight extrapolation of his results gives T,=0.19 
x10-'° sec. Because of the complexity of the de- 
excitation process, no reliable estimate is avail- 
able for T,; however, some orientation may be 
gained by consideration of a related process. An 
equation analogous to Eq. (3) may be used to cal- 
culate the capture rate for = hyperons produced 
by K” absorption in hydrogen. The observed cap- 
ture rate is obtained with T, =(0.10+ 0.25) x107" 
sec.’ To relate this result to the present situa- 
tion, we make the rough assumption that the mass 
dependence is not worse than that characteristic 
of energy loss to electrons by collision," i.e., 


D_ HH -10 
- = (Hp / HT, ~(0.2+ 0.5) x10" sec, 


Table II. Decay rate of £™ relative to the free 2. 
The =” was assumed to be bound by 0.5 Mev. 





Amplitude for 


x decay into Bound state of 2” 





S and P channels 1S 3S 
S#0; P=0 1.28 0.94 
S=0; P~0 0.94 1.04 








where Up (ug) is the reduced mass of the D~n-d 
(=~-p) systems. 

Under the above assumption, Eq. (3) gives P, 
= 0.327?!" for the most unfavorable case (K~-> 
parity even, = m bound in the ’S state). Using 
this value, we estimate that if the fraction of > 
productions leading to the bound state were 14%, 
6ts =n capture events should have been observed 
Since the entire distribution of events is com- 
pletely consistent with all = ’s being produced 
via reaction (2), and since, at most, one event 
can reasonably be interpreted as a hyperfragment 
production, it is probable that the branching ratio 
for reaction (1) is substantially less than 1%. 

We would like to thank Professor A. Pais and 
Professor George Snow for informative discus- 
sions during the course of this work. It is a 
pleasure to acknowledge the advice and encourage- 
ment of Professor Luis Alvarez and the assist- 
ance of the bubble-chamber and scanning staffs oj 
Lawrence Radiation Laboratory. 
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SPHERICAL GRAVITATIONAL WAVES* 


Ivor Robinson 
Department of Physics, University of North Carolina, Chapel Hill, North Carolina, 


A. Trautman 


Institute of Physics, Polish Academy of Science, Warsaw, Poland 


+ 


(Received March 7, 1960; revised manuscript received March 24, 1960) 


This note presents a class of solutions to Ein- 
stein’s gravitational equations for empty space. 
Some of the solutions appear to represent a very 
simple kind of spherical radiation. The metric 
considered has the form 


ds* = 2dpdo+(K - 2Hp - 2m/p)do? 
- pp-*{|dé + (8q/8n)do)? + [d+ (2q/08)do}*}, 


where m is a function of o only, p and gq are 
functions of o, &, 7, 


H=p"' ap/d0+ pe*p-*g/ tan - pq p~*/dé on, 


and K is the Gaussian curvature of the surface 
p=1, o=constant, 


K = p?(87/8&? + 87/87?) lnp. 


for this metric, the empty-space condition Rj, 
=0 reduces to 


8°q/8&? + &q/an? = 0 


(from which we can derive g=0 by a coordinate 
transformation), and 


0K /8&* + 8K / ar? = 4p-7(8/80 - 3H)m. 


If these equations are satisfied, the curvature 
tensor may be written as 


ence -2 -1 
Rig =P Dig tO gn tO Nip 


where Dijri, Wijk, and Njjp] are tensors of type 
Idegenerate, type III, and type II null, respect- 
ively. They are covariantly constant on any ray 
of constant o, é, 7. 

The solution is degenerate type I if m is non- 
zero and K is independent of é, 7. It is then re- 
ducible to m=1, p=coshut, g=0, where up isa 
teal or purely imaginary constant. If y is real 
and nonzero, this is Schwarzschild’s solution for 
amass p~*, 

If (€K/8¢)? + (aK /8n)? + 0 and the empty-space 
‘quations are satisfied, the metric is type II non- 
tull or type III, the condition for type III being 
m=0. As an example of these types, we mention 
m=constant, p= *?, qg=0. This contains a re- 


gion of incorrect signature, which can be elimi- 
nated, however, by the substitution ¢ =4¢~*. 

In the remaining case, where m vanishes and 
K is independent of £, 7, the solution is type II 
null or flat, the condition for flatness being A =0, 
where 


A =(8/8£ + i8/8n)[ p?(8/dE + 18/8n)H |. 


For K nonzero, |A/Kp| is invariant under all 
coordinate transformations which preserve the 
form of the metric. It might be described as the 
intensity of the gravitational field. 

This metric and the plane gravitational wave 
have the same local geometry, since their cur- 
vature tensors are algebraically indistinguishable. 
The wave fronts are subspaces of constant o. In 
the linear approximation, they are spheres ex- 
panding at the speed of light, the source having a 
velocity less than, equal to, or greater than this, 
according as K is positive, zero, or negative. 

The coordinates ~ and 7 may be chosen so that 
p=1+4(t?+7?) K(o). All components of the cur- 
vature tensor-—covariant, contravariant, and 
mixed—are then homogeneously linear ing. The 
contravariant metric tensor may be divided into 
a flat background metric, independent of g, and 
a residue which is homogeneously linear in q. 
There is thus an uncommonly close relation be- 
tween the rigorous solution and its linear approxi- 
mation. 

It is easy to construct nontrivial solutions which 
are periodic in o. In such a solution, there is no 
monotonic change which could be identified as 
loss of energy by the source. For positive K, 
however, there is at least one singularity on any 
wave front where the field does not vanish iden- 
tically; and P. G. Bergmann has suggested that 
this might conceivably represent a flow of matter 
which restores to the source the energy carried 
away by radiation. 

We can derive the plane-fronted wave,’ 


ds? = 2dp do - 2Hdo? - dé? - dr}, 
8H / dp = 8°H /dt? + 87H / 87? =0, 
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from this solution by specializing p as before, 
substituting p=2-?+r-"f, o=AB, E=r7E, N=r7h, 
q=)*9, where \ is constant, and taking the limit 
as \ tends to zero. There is then a singularity 
on every wave front where the homogeneity con- 
ditions 8°H/8é* = 88H /8n*> =0 are violated. 

We are much indebted to Professor P. G. Berg- 
mann for his comments on an earlier draft of 


this note. 
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POISSON BRACKETS BETWEEN LOCALLY DEFINED OBSERVABLES IN GENERAL RELATIVITY* 


Peter G. Bergmann and Arthur B. Komar 
Department of Physics, Syracuse University, Syracuse, New York 
(Received March 28, 1960) 


In this note we describe a new method for 
formulating local observables, and Poisson 
brackets between them, in Einstein’s theory of 
gravitation. The applicability of this method 
depends on the functional independence of the 
four scalars of the Weyl tensor but involves no 
global or topological assumptions. It combines 
the Hamiltonian approaches of Dirac’ with the 
construction of observables by Komar,’ and it 
leads to closed-form expressions both for the 
observables and the Poisson brackets. 

Our point of departure is the discovery that the 
four scalars of the Weyl tensor can be expressed 
in closed form in terms of the 12 canonical vari- 
ables g,,»,>””” as defined by Dirac. Using the 
notation, 


y= (g-v2) Of (1) 
mn 


mn 


to denote Dirac’s “invariant velocities,” we find 
for the components of the Weyl tensor the follow- 
ing expressions: 


~ re 
I a dd) ed af 


Pp _ J 
Cant inio! =" ise’ efi’ 


p,o im 
Cn Corto! hada iklm® (2) 
°Rikim denotes the three-dimensional curvature 
tensor, and the solidus signifies covariant dif- 
ferentiation with the help of three-dimensional 
Christoffel symbols. The Weyl scalars are 
quadratic and cubic expressions in these com- 
ponents. For instance, the first two (the quadra- 
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tic) scalars are 


1 iklm ikl kl 

= iktm© +4C€ 1° +4C, ; 

ik k 

Coa” ote” of 6) 
e/MS is Levi-Civita’s (three -dimensional) fully 
antisymmetric tensor, and indices are raised 
and lowered with the help of the three-dimen- 
sional metric e””, g,,». 

In order to retain flexibility, we shall assume 
that the four intrinsic coordinates to be used will 
be some four functions f?(A!, ..., A*) whose 
specification we may reserve; if desired they 
may be chosen so that the f coordinates are 
asymptotically Lorentzian for some specific 
Riemann-Einstein manifold. We shall, accord- 
ingly, introduce the four coordinate conditions 


A 


f’ -x" <0, (4) 


which determine the coordinate system uniquely. 
These four coordinate conditions, along with the 
constraints 


(5) 


form a system of eight second-class constraints, 
and hence lend themselves to the construction 

of field variables and functionals whose Poisson 
brackets with all the constraints (4), (5) vanish. 
Such variables are then observables in the techni- 
cal sense in which this expression is used now 
in general relativity. They are identical with the 
observables constructed in reference 2, but now 
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they appear as closed-form expressions in terms 
of Dirac’s 12 canonical field variables. 

As a preliminary we note that the remaining 
components of the metric tensor are uniquely 
jetermined by the coordinate conditions, and that 
they are closed-form expressions in terms of the 
canonical rrr We have 


(g° 01/2, (,x")=[F°,H,'1, 


ge" ow,x) =e ya). 8) 


for the consistency of the whole scheme it is of 
the essence that the Poisson brackets on the right 
contain only the delta function as a factor, and 

not, e.g., derivatives of the delta function. This 
condition is satisfied automatically, because the 
constraints (5) are generators of infinitesimal 
coordinate transformations, and the f?, being 
scalars, have transformation laws that consist 
solely of the transport term, 


p 


yy” - 3” ae (7) 


Since the right-hand side contains the 5x” only 
undifferentiated, no derivatives of the delta func- 
tion will occur on the right of Eqs. (6). 

To construct an observable we start from any 
functional of the canonical variables, A. We shall 
now add to A a linear combination of the con- 
straints (4), (5), so that the new functional A* 
commutes with all of them, 


A*- A+Sty -x )+€ H. +eHH la’. (8) 


Provided the constraints are satisfied, A* will 
equal in value the original A. The coefficients y, 


€ are determined as follows: 


Yp =[H,,4], 


00,1/2 


Os 
Y= [H,»A]-s Ys? 


_tk MF, eh, 
=-[f Al+ fi m ly,‘d a 


Pas 1/2. L 


Ale fee ae.) 
The inversion of the matrix of the Poisson brac- 
kets of the second-class constraints is thus 
accomplished locally and without difficulty. The 
Poisson brackets between any two observables 
A*, B* can now be calculated by means of the 
conventional Poisson brackets between canonical 
field variables. 

The Poisson bracket between A* and the Hamil- 
tonian vanishes, because the Hamiltonian is but 
a linear combination of the constraints (5), as 
was originally shown by Dirac.’ Nevertheless, 
A* is not a constant of the motion, but rather we 
have 

* - 

A fuller account of our results will be sub- 

mitted to the Physical Review. 
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ABSTRACTS 








In this section are printed the abstracts of Articles 
that have been forwarded to The American Institute of 
Physics for publication in THE PHYSICAL REVIEW. 
In quoting information obtained from this section be- 
fore the appearance of the corresponding Article, ref- 
erence should be made to “Physical Review (to be 
published)” rather than to this Journal. 


IRREVERSIBLE THERMODYNAMICS OF A 
NONLINEAR R-C SYSTEM. William Bernard, 
Research Division, Raytheon Company, Wal- 
tham, Massachusetts, and Herbert B. Callen, 
Department of Physics, University of Pennsyl- 
vania, Philadelphia, Pennsylvania (Received 
January 26, 1960). 


The formalism of nonlinear irreversible 
thermodynamics, as recently given by us, is 
applied to a circuit containing a capacitor and a 
nonlinear resistor. The solution is compared 
with those obtained by Lax and other investiga- 
tors. It is shown that the fluctuation-dissipation 
theorem is rigorous, and that no correction 


factors need be introduced for nonlinear systems. 


The dynamical behavior of the microscopic fluc - 
tuations, from which the macroscopic motion 
can be obtained, is also derived. Finally, a 
specific Markoffian model of a nonlinear R-C 
system, in strong interaction with a temperature 
reservoir, is shown to be consistent with the 
general analysis. 


IONIC CONDUCTIVITY OF HIGHLY IONIZED 
PLASMAS. M. Sakuntala, A. von Engel, and 
R. G. Fowler, Physics Department, University 
of Oklahoma, Norman, Oklahoma (Received 
March 30, 1959; revised manuscript received 
March 17, 1960). 


When a cloud of highly ionized gas ejected by 
a plasma shock tube is made to travel across a 
constant magnetic field, an electromotive force 
is produced in the plasma in a direction normal 
to both the field and the plasma path. Using two 
probes facing one another this electromotive 
force has been measured with an oscillograph. 
Its maximum value was found to be proportional 
to the field and the probe separation. By taking 
the maximum probe potential for different values 


434 





ee 


of the external resistance between the probes, 
the lowest value of the “resistivity of the plasma’ 
as measured by a current entering and leaving 
it was obtained. The resistivity has been shown 
to be independent of the magnetic field, the 
collecting area, the separation and surface state 
of the probes. All experiments were made in 
hydrogen at a gas pressure between 0.5 and 5 
mm Hg with a nearly critically damped current 
pulse of order 10* amperes lasting for about 
6-8 psec and fields< 2000 gauss. 

The plasma resistivity between the probes was 
found to be of the order of 1 ohm cm at a gas pres. 
sure of a few mm of Hg. This is about 100 times 
larger than the electronic resistivity of a fully 
ionized gas for currents circulating internally. 
The degree of ionization is thought to be high 
enough for interaction between charged particles 
to predominate. The measured values of the 
plasma resistivity agree with the results obtained 
from theory based on ionic conduction in the 
plasma which is here the necessary prerequisite 
for maintaining the electric neutrality of the 
moving plasma. 

From the measured probe voltage the flow ve- 
locity of the plasma was derived. Its variation 
with gas pressure agrees with shock-wave theory. 


MEANING AND INTERPRETATION OF ACOUS- 
TIC MOMENTUM AND ACOUSTIC RADIATION 
STRESS. E. J. Post, Bell Telephone Labora- 
tories, Murray Hill, New Jersey (Received 
April 3, 1958; revised manuscript received 
May 20, 1959). 


This paper surveys viewpoints and available 
procedures of analysis of acoustic radiation 
stress and the concept of acoustic momentum. 

An alternative approach is suggested based on 
Eckart’s principle and Poynting’s condition of 
separate observability of matter and elastic wave 
motion. The results are consistent with those 
obtained by means of the adiabatic principle but 
show discrepancies with the so-called Eulerian 
approach for transverse waves. 

Radiation stress and radiation momentum are 
typical field concepts needed to complement the 
continuum description. In a sense they serve to 
establish consistency with the particle approach, 
thus filling a need which the conventional theory 
of elasticity does not satisfy. 

The energy-momentum tensor of elastic wave 
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motion provides a unique example of a physically 
meaningful space-time covariance in the frame- 
work of Galilean kinematics. 


OPTICAL PROPERTIES OF COPPER. S. Roberts, 
General Electric Research Laboratory, Schenec- 
tady, New York (Received December 21, 1959). 


Optical properties of solid copper have been 
measured in a wavelength range from 0.365 to 
2.5 microns and at temperatures 90°, 300°, and 
500°K. Annealed and electropolished specimens 
were used. Measurements were done ina 
vacuum of about 21075 mm (Hg) and it was ob- 
served that the oxide film which normally forms 
on copper upon exposure to air could be removed 
by heating in vacuum to a temperature of 500°K. 
In the wavelength range below 0.6 micron the 
data confirm the well-known facts relating to a 
threshold for interband electronic transitions. 

At longer wavelengths the optical properties are 
determined almost entirely by free electrons. 
Deviations from simple theory are partly ex- 
plained by the anomalous skin effect and it is 
concluded that the electronic collisions at the 
metal surface are diffuse. However, it is shown 
that the anomalous skin effect is not sufficient 

to explain the observed deviations and that a 

more complete interpretation ought to consider 
the nonspherical nature of the Fermi surface and 
variation of the relaxation time over this surface. 


SPECIFIC HEAT OF DILUTE ALLOYS. W. Mar- 
shall, Theoretical Physics Division, Atomic 
Energy Research Establishment, Harwell, 
England (Received January 4, 1960; revised 
manuscript received March 1, 1960). 


Recently Zimmerman has observed that the 
addition of Mn to Cu produces a large contribu - 
tion to the specific heat which, at low tempera- 
tures, is linear in temperature and independent 
of Mn concentration. It is shown that: (1) this 
remarkable result can be explained in terms of 
the well-known Ruderman-Kittel-Yosida spin- 
Spin coupling via conduction electrons; (2) the 
specific heat results of Beck et al. on FeV and 
FeCr alloys are probably of essentially the same 
origin as those of Zimmerman on CuMn; and (3) 
there are serious objections to the mechanism 
of antiferromagnetism postulated by Overhauser 





and used by him to explain the specific heat re- 
sults. In contrast to the Overhauser theory, no 
new concepts are involved and it is suggested 

that the large specific heat comes from a small 
fraction of Mn spins which, because of the ran- _ 
dom nature of the alloy, happen to be in small 
effective fields and therefore not strictly aligned. 
The theory depends on two plausible assumptions 
which have not, at present, been proved rigorous- 
ly valid. 


NUCLEAR SPIN-LATTICE RELAXATION 
CAUSED BY PARAMAGNETIC IMPURITIES. 

W. E. Blumberg, Department of Physics, Uni- 
versity of California, Berkeley, California (Re- 
ceived December 24, 1959). 


The behavior of the nuclear magnetization has 
been examined for the transient region immedi- 
ately following saturation of the nuclear spin 
system in the case in which the nuclear relaxa- 
tion is limited by spin diffusion to paramagnetic 
impurities. Relaxation effects in the absence of 
diffusion and the presence of rapid diffusion are 
discussed. An experiment is reported which 
substantiates the calculations for the case of 
diffusion-limited relaxation and, in addition, 
allows a determination of the correlation time 
of the magnetic moment of the paramagnetic ion 
and the coefficient of nuclear spin diffusion. 


g FACTOR AND DONOR SPIN-LATTICE RE- 
LAXATION FOR ELECTRONS IN GERMANIUM 
AND SILICON. Laura M. Roth, Lincoln Labora- 
tory, Massachusetts Institute of Technology, 
Lexington, Massachusetts (Received January 
29, 1960). 


The g factors of electrons in Ge and Si are 
calculated on the basis of the effective-mass 
approximation. The results are consistent with 
experimental spin resonance data. The effect is 
predicted to be anisotropic, with g i less than 
g,- This anisotropy introduces a strong inter - 
action between the electron spin and shear 
waves, for the singlet donor ground state. This 
interaction can account in order of magnitude 
for the observed spin-lattice time 7, for donor 
electrons in Si at low temperatures, including 
both a one-phonon process and a two-phonon 
Raman-type process. The temperature and mag- 
netic field dependence for the two processes 
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are predicted to be 7,~* « TH* and 7,™* « T"H’, 
respectively. The temperature dependence 
agrees with experiment; however, there are 
discrepancies in the magnetic field dependence. 
T, is predicted to be anisotropic, an effect that 
has not as yet been observed. 


IRREVERSIBILITY IN THE SUPERCONDUCTING 
TRANSITION OF LEAD. R. W. Shaw and D. E. 
Mapother, Physics Department, University of 
Illinois, Urbana, Illinois (Received January 20, 
1960). 


The nature and origins of the hysteresis in the 
magnetic transition of pure superconducting lead 
have been studied in some detail. The effect can 
be introduced in nonhysteretic specimens by 
plastic deformation at liquid helium tempera- 
tures. The hysteresis becomes apparent some- 
what below T. and its width increases monoto- 
nically with decreasing temperature. Similar 
effects can be caused by dilute additions of Ca 
to Pb. Strain-induced hysteresis anneals out 
near 300°K but temperatures approaching the 
melting point appear necessary to remove 
impurity- induced effects. Isothermal resistive 
measurements show a small fraction of the 
superconducting phase persisting to fields 
several hundred gauss above H,. The residual 
superconductivity at high fields is increased by 
plastic deformation and, in general, seems 
closely related to the hysteresis effect. The ob- 
servations suggest the existence of a connected 
network which pervades the entire specimen 
volume and consists of very small filaments 
having a critical field exceeding the reversible 
critical field of bulk lead. The filaments are 
believed to be associated with defects in the 
crystalline lattice. 


PARAMAGNETIC RESONANCE AND OPTICAL 
SPECTRA OF DIVALENT IRON IN CUBIC 
FIELDS. I. THEORY. W. Low and M. Weger, 
Department of Physics, Hebrew University, 
Jerusalem, Israel (Received January 22, 1959; 
revised manuscript received January 20, 1960). 


The energy level splittings of the ground state 
of the d® configuration in cubic and axial fields 
are given. The Zeeman splittings of the various 
levels are calculated for weak and strong mag- 
netic fields. In the case of tetrahedral symmetry 
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the effect of the perturbations of the odd-parity 
configurations d°p and d*f on the ground state is 
estimated. 


PARAMAGNETIC RESONANCE AND OPTICAL 
SPECTRA OF DIVALENT IRON IN CUBIC 
FIELDS. II. EXPERIMENTAL RESULTS. 

W. Low and M. Weger, Department of Physics, 
Hebrew University, Jerusalem, Israel (Re- 
ceived January 22, 1959; revised manuscript 
received January 20, 1960). 


The paramagnetic resonance absorption of 
Fe** in MgO is observed at g=3.428 and 6.86. 
The optical absorption line is found at 10000 
cm~*. The paramagnetic resonance spectrum 
indicates considerable covalent bonding. The 
origin of the line at 6.86 is discussed. 

In tetrahedral ZnS a paramagnetic line is 
found at g=2.25 and optical absorption at 3 yu 
and 0.7 ». Possible explanations of this spec- 
trum are discussed. 

A short discussion of the optical absorption 
spectra of trivalent iron in MgO is presented. 


INFLUENCE OF F CENTERS ON THE LATTICE 
THERMAL CONDUCTIVITY IN LiF. Robert 0. 
Pohl, Laboratory of Atomic and Solid-State 
Physics, Cornell University, Ithaca, New York 
(Received January 25, 1960). 


The influence of photochemically produced 
F centers in LiF on the thermal conductivity is 
investigated, the density of the F centers being 
determined optically. The F centers decrease 
the thermal conductivity appreciably at low 
temperatures. Additive coloration of KCl has 
a similar effect. The experimental results are 
quite different from the results obtained in the 
case of the isotope effect in Ge. Callaway’s 
theory for the lattice thermal conductivity in the 
presence of point defects cannot explain the 
present observations satisfactorily, although it 
is superior to that of Klemens. Two explana- 
tions are proposed: (1) The F centers are not 
randomly distributed (but the presence of clus- 
ters of F centers can be excluded in the experi- 
ments). (2) The long-range strain field around 
the F center does not scatter like a point defect. 
In support of Callaway’s theory it is shown that 
his model can explain an experimental observa- 








arity 
ate is 


ICAL 


/sics, 


ion 
ted, 


[TICE 
rt O. 


York 








VoLUME 4, NUMBER 8 





PHYSICAL REVIEW LETTERS 


APRIL 15, 1960 





tion (first reported by Toxen) about the influence 
of point defects on the thermal conductivity. 


THERMAL EQUILIBRIUM OF COLOR CENTERS 
IN DOPED KCl CRYSTALS. P. Camagni,* S. 
Ceresara,* and G. Chiarotti, Istituto di Fisica 
dell’Universita, Pavia, Italy (Received March 2, 


1959; revised manuscript received June 15, 1959). 


In additively colored crystals of KCl:Sr a 
thermal reaction between F and Z, centers is 
observed. The reaction is reversible and ap- 
proaches equilibrium in a few hours, in the 
range of temperatures between 70°C and 170°C. 
The equilibrium value of the ratio between the 
concentrations of F and Z, centers, (xp/xz,), 
does not depend upon the total amount of colora- 
tion nor upon the past history of the crystal, and 
for a given impurity content is a function only 
of temperature. For crystals containing 2.5 x10~* 
molar fractions of Sr, the experimental points 
fit the equation 


x_/x_ =0.8x10~*exp(-0.28/kT), kT in ev. 


FZ, 
A comparison of the empirical equation with 
those that can be derived through the application 
of the mass-action law to some specific models 
of the F-+ Z, conversion is discussed. 


"Now at Centro Nazionale di Studi Nucleari, Ispra, 
Italy. 


MAGNETISM OF EUROPIUM GARNET. W. P. 
Wolf, Clarendon Laboratory, Oxford, England, 
and J. H. Van Vleck, Harvard University, Cam- 
bridge, Massachusetts (Received January 18, 
1960). 


The theoretical expressions for the magnetic 
moment of a trivalent europium ion in a molecu- 
lar field arising from exchange are applied to 
Pauthenet’s measurements on europium iron 
garnet. It is a good approximation to assume 
that the exchange interaction stems entirely 
from the coupling with the iron atoms, which 
greatly simplifies the theory since the molecu- 
lar field on the europium is then an impressed 
one and does not have to be determined self- 
consistently. The calculated variation of the 
magnetization with temperature is in excellent 
accord with experiment. The magnitude of the 
exchange interaction is compared with that in 









the other rare earth iron garnets; it is almost 
exactly the same as in gadolinium iron garnet. 


MAGNETIC SUSCEPTIBILITY OF METALLIC 
EUROPIUM. R. M. Bozorth, Bell Telephone 
Laboratories, Murray Hill, New Jersey, and 

J. H. Van Vleck, Harvard University, Cambridge, 
Massachusetts (Received January 21, 1960). 


The susceptibility of europium is measured 
from 1.3 to 300°K in fields up to 12 000 oersteds. 
This metal is found not to be ferromagnetic, but 
to have at low temperatures a very high para- 
magnetic susceptibility, about forty times higher 
than for the free ion or hydrated salts. Varia- 
tion of susceptibility with field strength is ob- 
served below 100°K. The saturation curvature 
at low temperatures is very large, and practi- 
cally independent of temperature. The suscepti- 
bility at high temperatures is consistent with a 
divalent model. The magnetic behavior at low 
temperatures is hard to interpret on the basis 
of a divalent ion and can be more readily ex- 
plained with a trivalent one for which the theory 
is developed. Metallurgical evidence, however, 
indicates that metallic europium is divalent 
even at low temperatures. On the other hand, 
it is generally believed that europium is triva- 
lent in Eulr,, but if one uses the conventional 
molecular field theory the reported ferromag- 
netism of Eulr, cannot be ascribed to a Eu*** 
with a reasonable value of the exchange field 
consistent with that in Gdlr,. 


ELASTIC CONSTANTS OF PALLADIUM FROM 
4.2 TO 300°K. J. A. Rayne, Westinghouse Re- 

search Laboratories, Pittsburgh, Pennsylvania 
(Received January 20, 1960). 


Elastic-constant measurements have been 
made on a single crystal of palladium in the 
temperature range 4.2-300°K. Extrapolation of 
the data to absolute zero gives 


Cy, = (2.341 + 0.027) x10"? dyne cm~, 
Cyp = (1.761 + 0.027) x10" dyne cm~, 
Cag = (0.712 + 0.003) x10” dyne cm~. 


The corresponding value of Debye temperature 
is 6,=275+8°K, which compares well with the 
calorimetric figure of 6,=274+3°K. Both shear 
constants show an anomalous temperature de- 
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pendence. This dependence can be correlated 
with the temperature variation of the contribu- 
tion to C=c,, and C’=3(c,,-¢,,), resulting from 
the presence of holes in the d band of palladium. 


SPIN-LATTICE RELAXATION OF SHALLOW 
DONOR STATES IN Ge AND Si THROUGH A 


DIRECT PHONON PROCESS. Hiroshi Hasegawa,~ 


Carnegie Institute of Technology, Pittsburgh, 
Pennsylvania (Received January 22, 1960). 


The many-valley character of the conduction 
band edge of germanium and silicon causes an 
anisotropy of the g shift and of the deformation 
potential for the conduction electrons. It is 
shown that the combination of these two effects 
provides a mechanism for spin-lattice relaxa - 
tions of the donor spins in germanium and sili- 
con that yields 1/7, proportional to the tem- 
perature T and to the fourth power of the static 
magnetic field H. Using known data about the 
deformation potential constant, the g shift, the 
energy of the intervalley splitting, and the 
elastic constants, the magnitude of T, is found 
to be approximately 2 x10~° sec for phosphorus 
donors in germanium, and 1 x10* sec for phos- 
phorus donors in silicon. These values refer 
to T =1.25°K, H =3000 gauss, with the field 
applied along the |111| axis. Our mechanism 
fails to give a finite 7, for donors in silicon, 
when the field is applied along the [100] axis. 


*On leave of absence from the Department of Physics, 
the University of Tokyo, Tokyo, Japan. 


MAGNETIC PROPERTIES OF THULIUM METAL. 
D. D. Davis and R. M. Bozorth, Bell Telephone 
Laboratories, Murray Hill, New Jersey (Re- 
ceived January 21, 1960). 


Measurements made from 1.3 to 300°K in fields 
up to 12000 oe, show that thulium becomes anti- 
ferromagnetic on cooling to 60°K and is ferro- 
magnetic below 22°K. In agreement with earlier 
work the susceptibility between 70 and 300°K 
obeys the Curie-Weiss law and the effective 
Bohr magneton number is 7.6, equal to the 
theoretical value for the trivalent ion core hav- 
ing two holes in the 4f shell. The highest ob- 
served ferromagnetic moment is only 0.5 Bohr 
magneton whereas the theoretical value is g@/J =7. 
Large hysteresis and crystal anisotropy are ob- 
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served at the lowest temperatures. It is sug- 
gested that the anisotropy is so large that the 
magnetization cannot be deviated appreciably 
from the hexagonal axis even in the higher fields 
(70000 oe) used by Henry; he observed 3.4 Bohr 
magnetons whereas the highest expected for poly. 
crystalline material on this assumption is 3.5, 


DETERMINATION OF THE EFFECTIVE SCAT. 
TERING MECHANISM PARAMETER OF ELEC. 
TRON TRANSPORT THEORY. Theodore C. 
Harman, Lincoln Laboratory, Massachusetts 
Institute of Technology, Lexington, Massachu- 
setts (Received January 25, 1960). 


For Corbino disk geometry, the variation of 
thermoelectric power, a, with magnetic field 
strength, H, is a strong function of the scatter- 
ing and a weak function of the Fermi level. The 
zero magnetic field thermoelectric power isa 
strong function of the Fermi level and a weak 
function of the scattering. Thus, measurements 
of Aa/a(0) versus H and a(0) at a particular 
impurity concentration and temperature should 
yield a scattering parameter and a Fermi level 
separately. The transport equation for Aa/a(() 
has been derived for the case of a simple 
spherically symmetric energy band, arbitrary 
scattering, exact statistics, and Corbino disk 
geometry. The thermoelectric power increase 
55% as the magnetic field increased to 22 kgauss 
for mercuric selenide. Although there is some 
evidence which indicates that the simple band 
form does not occur in HgSe, a preliminary 
comparison of theory and experiment showed 
satisfactory agreement. 


INTERNAL FIELD EMISSION AT NARROW 
p-n JUNCTIONS IN INDIUM ANTIMONIDE. 

A. G. Chynoweth and R. A. Logan, Bell Tele- 
phone Laboratories, Murray Hill, New Jersey 
(Received January 25, 1960). 





An experimental study has been made of the 
field and temperature dependence of internal 
field emission in narrow p-n junctions in indium 
antimonide. Relatively good agreement, both 
qualitative and quantitative, is obtained between 
the experimental results and the usual expres- 
sion for the barrier transparency. From studies 
of Esaki characteristics at low temperatures 
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and from the observed temperature dependence 

of the tunnelling current, it is confirmed that 

the tunnelling transitions do not involve phonons. 
Also, it is shown that the temperature dependence 
of the barrier transparency is determined by 

that of the energy gap at k=0. 


ELECTRON PARAMAGNETIC RESONANCE OF 
Fe* IN TiO, (RUTILE). David L. Carter and 
Akira Okaya, Department of Physics, Columbia 
University, New York, New York (Received 
January 20, 1960). 


Paramagnetic resonance has been observed 
from 2 kMc/sec to 110 kMc/sec in Fe**-doped 
single crystals of rutile at 78°K, 4.2°K, and 
1.4°K. The rather large zero-field splittings 
measured between the three doublets are 43.3 
+0.1kMc/sec and 81.3+0.1 kMc/sec. The de- 
rived constants for a spin Hamiltonian are 
D=20.35+ 0.1 kMc/sec, E =2.21+0.07 kMc/sec, 
a=1.1+0.2 kMc/sec, F = -0.5+0.3 kMc/sec, 
g=2.000+ 0.005. The average linewidth is 35 
Mc/sec for 2x10" spins per cm', and the aver- 
age spin-lattice relaxation time T, at T=1.4°K 
is4 msec. T, still has a very slow inverse tem- 
perature dependence even at T=78°K. An in- 
crease in sensitivity over metal-walled cavity 
spectrometers was achieved by using the pieces 
of Fe-doped rutile as the microwave resonators. 


RESONANT ELECTRON CAPTURE AND STRIP- 
PING IN MODERATELY LARGE-ANGLE ATOMIC 
COLLISIONS. F. P. Ziemba,”* G. J. Lockwood, 
G.H. Morgan, and E. Everhart, Physics De- 
partment, University of Connecticut, Storrs, 
Connecticut (Received January 21, 1960). 


Differential scattering of ions by atoms in the 
energy range of 1 kev to 200 kev has been studied. 
The incident ion, after a single collision which 
is hard enough to result in a 5° deflection, was 
analyzed to determine whether it had captured or 
lost electrons. The angle 5° was chosen as a 
typical moderately large angle and held fixed as 
the energy of the incident ion was varied. When 
the electron capture probability is plotted vs 
energy, one or more peaks are observed. For 
the symmetrical case of He* on He, seven peaks 
are clearly outlined. Four peaks appear in the 





H* on He combination and three with H* on H,. 
Single or double peaks are found in other cases 
studied which include H* on N,, O,, air, Ne, Ar, 
and Kr; H,* on H, and He; H,* on H, and He; 

He* on Ne, Ar, and Kr; N* on Ar; N,* on N,; 
Ne* on Ne and Ar; Ne** on Ne and Ar; Ar* on 
Ar; and Kr* on Kr. For each case the probabil- 
ities for electron capture, scattering without 
change of charge, and various degrees of elec- 
tron stripping are plotted vs energy. In those 
cases in which the electron capture probability 
curve has more than two peaks, these peaks are 
nearly evenly spaced when the probabilities are 
plotted vs the time of the interaction. This in- 
dicates an electron exchange effect whose period 
is of the order of 10°*® second. In cases where 
there are many electrons involved in the colliding 
atoms the phenomenon is more complicated, but 
vestiges of this resonant exchange are sometimes 
observed. 


*Now at Hughes Aircraft Company, Culver City, 
California. 


APPROACH TO EQUILIBRIUM IN QUANTAL 
SYSTEMS: MAGNETIC RESONANCE. A Sher* 
and H. Primakoff, Washington University, 

St. Louis, Missouri (Received January 28, 1960). 


The present paper presents a derivation of the 
“master” or Boltzmann “gain-loss” equation 
from the Schrédinger equation, i.e., a deriva- 
tion of the equation for the evolution in time of 
the probabilities of finding a physical system 
in its various states from the equation for the 
corresponding probability amplitudes. The 
“master” equation is derived for an, in effect 
completely self enclosed, “supersystem,” 
[A+B], consisting of a “system of interest,” 
[A], and a “surroundings,” [B], in relatively 
weak mutual interaction. A discussion is given 
of the range of validity of the “master” equation 
for [A +B] and it is shown that the random phase 
assumption is required for the state vector of 
[A +B] at the initial time only. The normally 
microcanonical character of the equilibrium 
statistical configuration of [A +B] is demon- 
strated and a treatment is given of exceptional, 
“extremely quantal-coherent,” initial statistical 
distributions of [A +B] which may evolve away 
from equilibrium. Derivations are also pre- 
sented of the “master” equation for [A | and of 
the “master” equation for an individual particle 
or quasi-particle [q], within [A]; a discussion 
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of the range of validity of these “master” equa- 
tions is given and the normally canonical char - 
acter of the equilibrium statistical configuration 
of [A] is deduced. General solutions of the 
“master” equations for [A+B], [A], and [q] are 
worked out and the relation between the princi- 
ples of “microscopic reversibility” and “detailed 
balance” and the nonoscillatory character of the 
approach to equilibrium are exhibited. A 
theorem is presented regarding the time varia- 
tion of the entropy of [A]. 

As illustrations of the general methods de- 
veloped, two important processes in magnetic 
resonance —the time variation of the longitudinal 
magnetization, ¢ u),, and the time variation of 
the transverse magnetization, (hss —are dis- 
cussed in some detail. It is shown that the vari- 
ation of (u), with ¢ and of (u’), with ¢ for a 
“nonrigid” lattice can be described by means of 
the “master” equation for an individual spin [q] 
and several special cases are discussed on the 
basis of the evaluation of the appropriate transi- 
tion probabilities; a comparison with the “spin- 
temperature” procedure is also appended. On 
the other hand, it is demonstrated that for a 
“rigid” lattice no description of the variation of 
(Ku), with ¢ can be given on the basis of a 
“master” equation; in this case, quantal coher - 
ence effects neglected in the derivation of the 
“master” equation from the Schrédinger equation 
are vital and (..’); must be evaluated by a rig- 
orous calculation of Trace{[appropriate time- 
dependent density matrix ]’}. 

*Now at Centre d’Etudes Nucléaires de Saclay, 
Gif-sur-Yvette, Seine et Oise, France. 


RETARDATION AND K X RAY RELATIVE IN- 
TENSITIES. G. R. Taylor and W. B. Payne, 
Virginia Polytechnic Institute, Blacksburg, 
Virginia (Received December 7, 1959). 


The matrix elements for electric dipole transi- 
tions with retardation have been numerically eva- 
luated for transitions between the 1S and higher 
discrete states of a single Dirac electron ina 
Coulomb field. The results are used to calculate 
the relative intensities of the principal K x-ray 
lines in elements of high atomic number. Com- 
parison of the retarded and nonretarded matrix 
elements confirms the earlier result for lead 
that the effect of retardation is significant in 
heavy elements only when the total angular mo- 
mentum of the electron changes. Comparison 
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with the experiments of Beckman indicate that 
this conclusion is correct, and that in the Ka, 
to Ka, intensity ratio the effect of retardation 
is more significant than screening effects. 





@ CORRECTION IN THE SPECTRA OF THE 
IRON GROUP. Giulio Racah and Yehuda Shadmj 
The Hebrew University, Jerusalem, Israel (Re- 
ceived January 29, 1960). 


A semiempirical correction term, proportion 
to the “seniority operator” Q, is introduced int 
the Slater formulas for the configurations 3d%4: 
of the second spectra of V, Cr, and Fe. The 
agreement with the experimental data, obtained 
by least squares, is improved by about 20%. 
The physical meaning of this correction is dis- 
cussed. 


EXPERIMENTAL VERIFICATION OF THE 
“INCOHERENT SCATTERING” THEORY FOR 
THE TRANSPORT OF RESONANCE RADIATIOK 
A. V. Phelps and A. O. McCoubrey, * Westing- 
house Research Laboratories, Pittsburgh, 
Pennsylvania (Received January 20, 1960). 


The predictions of the incoherent scattering 
theory of the transport of resonance radiation 
developed by Holstein and by Biberman are 
shown to agree with laboratory measurements 
of the decay constants for the intensity of the 
mercury 2537A resonance radiation following 
a period of optical excitation. Also, the relative 
importance of the diffusion of resonance atoms 
and the escape of resonance radiation as mech- 
anisms for the destruction of mercury atoms in 
the *P, or resonance state are determined from 
measurements of decay constants as a function 
of mercury density and cell radius. The ex- 
perimental results show that diffusion of the 
resonance atoms is negligible and that the pre- 
dictions of imprisonment theory are confirmed 
to within 15%. The experiment sets an upper 
limit to the diffusion coefficient at unit gas 
density for atoms in the *P, state of 5x10" cm 
sec™ at 340°K, which is consistent with a value 
of 4x10'® cm™ sec™ predicted using the fre- 
quency of excitation transfer collisions calcu- 
lated by Holstein. The success of the theories 
based on completely incoherent scattering of 
resonance radiation points to the desirability of 
including this feature in the treatment of those 
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astrophysical problems in which the spectral 
line shape is determined by Doppler broadening 
or by collision broadening. 


*Now at the National Company, Incorporated, Mel- 
rose, Massachusetts. 


ANALYSIS OF THE TWO-MODE-OF- FISSION 
HYPOTHESIS. George P. Ford, Los Alamos 
Scientific Laboratory, University of California, 
Los Alamos, New Mexico (Received March 16, 
1959; revised manuscript received March 3, 
1960). 


Mass yields for gamma fission of U***, proton 
fission of Th**?, alpha fission of U***, and deu- 
teron fission of normal uranium have been ex- 
amined in terms of vector spaces for consistency 
with the hypothesis that for each target and pro- 
jectile there are two and only two modes of fis- 
sion. The first three cases are on the whole 
consistent with the hypothesis. The fourth is not 
consistent with the hypothesis; the measured 
yields include some yields from neutron fission. 
Some consequences of the hypothesis with the 
added requirement of non-negative yields and 
non-negative coordinates are derived. 


RADIOCHEMICAL STUDIES OF THE (), pn) 
REACTION IN COMPLEX NUCLEI IN THE 80- 
450 Mev RANGE. Herbert P. Yule* and Anthony 
Turkevich, Enrico Fermi Institute for Nuclear 
Studies, University of Chicago, Chicago, IIli- 
nois (Received January 22, 1960). 


Excitation functions for the (p,pn) reaction have 
been determined in the energy range 82 Mev to 
426 Mev for the target nuclei F’®, Cu®*, and Au’”. 
The absolute values are based on the excitation 
function for the reaction C*(p,pn)C™“. These 
excitation functions exhibit a general decrease 
with energy and the cross sections lie between 
122 mb for Au’®* at 82 Mev and 23 mb for F"® at 
426 Mev. An excitation function for the reaction 
Al""(p, 3pn)Na** is also presented. The results 
of Monte Carlo nuclear cascade calculations 
have been used to predict these (p, pn) excitation 
functions. The theoretical results are compared 
with experimental results of this and other re- 
ports. Agreement is obtained for the F’*(p, pn) F*® 
excitation function and an extrapolation of the 
Monte Carlo results. The theoretical excitation 
functions are about one-half the experimental 





results for the Cu®(p,pn)Cu™ reaction and about 
one-third the experimental results for Au’”(p, 
pn)Au'**, The effects on the Monte Carlo calcula- 
tions of a variation in the radius parameter are 
examined. 


"Present address: California Research Corporation, 
La Habra, California. 


STRIPPING MECHANISM FOR REACTIONS WITH 
SMALL Q VALUE: THE REACTION Li’(d, p)Li®. 
J. P. F. Sellschop, Nuclear Physics Research 
Unit, University of the Witwatersrand, Johannes- 
burg, South Africa (Received December 8, 1959). 


A number of angular distributions of protons 
from the reaction Li’(d, p)Li® have been meas- 
ured for a range of incident deuteron energies 
below 2.5 Mev. These agree remarkably well 
with a simple form of Butler-Born stripping 
theory, uncorrected for Coulomb and nuclear 
effects. A description is given for this unusual 
agreement in terms of the small Q value, -0.188 
Mev, for the reaction. 

A resonance in the proton yield is found at an 
incident deuteron energy of 1.4 Mev which has 
not been observed in measurements of the 8 
yield from this reaction. Angular distributions 
measured on and around the resonance show no 
influence of this on the unusually good stripping 
patterns. 


GAMMA RAYS FROM DEUTERON STRIPPING 
REACTIONS. G. R. Satchler, Oak Ridge National 
Laboratory, Oak Ridge, Tennessee, and 

W. Tobocman, The Rice Institute, Houston, 
Texas (Received January 19, 1960). 


The distorted-wave Born approximation is 
used to calculate the p-) angular correlation 
from several deuteron-stripping (d, p) reactions. 
One /=2 and four /=1 captures are considered. 
Optical potentials with rounded edges are used 
to distort the wave functions. In some cases 
the correlation is considerably changed from 
the pattern predicted by the plane wave Born 
approximation, and the distortion effects are 
strongly dependent on the direction of the emitted 
proton, and on the type of distortion assumed. 
We include a general discussion of the theory of 
the (d, py) correlation. 
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STUDIES OF LOW-LYING LEVELS OF EVEN- 


EVEN NUCLEI WITH (d, p) AND (d, t) REACTIONS. 


Bernard L. Cohen and Robert E. Price, Radia- 
tion Laboratory, University of Pittsburgh, Pitts- 
burgh, Pennsylvania (Received December 3, 
1959). 


Low-lying states of even-even nuclei in the 
vibrational region were studied by exciting them 
with (d, p) and (d,t) reactions. The relative 
cross sections for exciting ground states (G) 
bear little relationship to whether they are al- 
lowed or forbidden by the simple Mayer -Jensen 
configurations, so that configuration mixing is 
generally large. Some of the details of this 
mixing are obtained. The allowed portions of 
these cross sections are generally quite close 
to the single-particle values. First and second 
excited states are much more strongly exciting 
than expected theoretically. In Pd’™, Pd*, 
Pt'*, and Pt’® searches for the triplet in second 
excited states indicate that its total spacing must 
be less than 80 kev. New states were found in 
the triplet region of Cd*"* and Cd"? bringing the 
total number of known states in these to 5 and 4, 
respectively; each includes two 0+ and two 2* 
states. Higher excited states were studied and 
in almost all cases they occur below the expected 
position of the third member of the vibrational 
band; this gives evidence on the size of the en- 
ergy gap. The location of all 0* levels up to 3 
Mev is determined for several nuclei. In two 
cases the 3* states required by the Davidow- 
Filipov theory are found to be 0*, and in other 
cases they are not found at all. A number of 
previously unknown levels are catalogued. 


PARAMAGNETIC AND OPTICAL SPECTRA OF 
YTTERBIUM IN THE CUBIC FIELD OF CAL- 
CIUM FLUORIDE. W. Low, Department of 
Physics, The Hebrew University, Jerusalem, 
Israel (Received January 28, 1960). 


The paramagnetic resonance spectrum of Yb** 
in CaF, was observed at 20°K and 3 cm wave- 
length. The spectrum is described by a cubic 
spin Hamiltonian = g8H -$+AS8-I with g£=3.426 
+0.001, S=1/2, A’ =(886.5+1.5) x10“ cm™, 
A‘"$ = (243.2+0.4)x10™ cm™, F'=1/2, P%=5/2. 
The ratio of magnetic moments is '/y'™ 
=1.374,+ 0.005. The optical spectrum shows 
lines of 9774 A, 9770 A, 9763 A, and more dif- 
fuse and unresolved bands at 9080 and 12730 A. 
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The paramagnetic spectrum is explained as 
arising from the [, doublet. The other levels 
are removed by at least a few cm™~ leading to 
an isotropic g value of 3A or 24/7 for the lowest 
I’, level. 


RESPONSE OF Nal(T1) TO ENERGETIC HEAVy 
IONS. E. Newman and F. E. Steigert, Yale Uni- 
versity, New Haven, Connecticut (Received 
January 22, 1960). 


The light output of Nal(Tl) was measured asa 
function of energy for ions of He*, B*°®, B™, Cc”. 
N*, O'*, F’®, and Ne” from approximately 1.0 
to 10.0 Mev/nucleon. The particle energy was 
varied with absorbing foils and the degraded 
beam analyzed in a magnetic spectrometer. 
Light output becomes linear with energy above 
approximately 6 Mev/amu. There appears to be 
an odd-even incident charge effect in the satura- 
tion value of the differential efficiency of fluores- 
cence, dL/dE. 


NUCLEAR SPIN AND HYPERFINE INTERACTION 
OF In™°”_ W. J. Childs and L. S. Goodman, 
Argonne National Laboratory, Lemont, Illinois 
(Received January 25, 1960). 


The hyperfine structure of the 1.7-hr, 393- 
kev metastable state of In*** has been studied 
by use of the atomic-beam magnetic -resonance 
technique. The nuclear spin is found to be 1/2 
and the magnetic dipole moment to be -0.21050 
+ 0.00002 nm, subject to a possible hyperfine 
anomaly. The hyperfine separation in the atomic 
P,,. state is measured to be 781.084 + 0.010 Mc/ 
sec. 


ANGULAR DISTRIBUTION OF FRAGMENTS 
FROM FISSION OF Au’*’ WITH CARBON IONS. 
Glen E. Gordon, Almon E. Larsh, and Torbjgrn 
Sikkeland, Lawrence Radiation Laboratory, 
University of California, Berkeley, California 
(Received January 26, 1960). 


The kinetic energy and angular distribution of 
fragments from fission of Au'®’ with 123- and 
93-Mev carbon ions have been determined by 
observation of the fragments in gas scintillation 
and solid-state detectors. Between 20 deg and 
160 deg in the center-of-mass system, both 
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angular distributions lie slightly above a 1/sin@ 
curve, falling below it beyond those angles. The 
anisotropies [0(0°)/o(90°)] are 4.7 and 3.8 for 
123- and 93-Mev carbon ions, respectively. 

The most probable fragment kinetic energies 

in the center-of-mass system are 73+3 Mev 

and 71+3 Mev. 


THEORETICAL PREDICTIONS FOR THE SPEC- 
TRA OF THE ODD-MASS XENON AND TELLU- 
RIUM ISOTOPES. Norman K. Glendenning, 
Lawrence Radiation Laboratory, University of 
California, Berkeley, California (Received De- 
cember 18, 1959). 


The spectra of the isotopes Te’**» 125, Xe!?7» 129 
are calculated by assuming that the odd neutron, 
having available the 3s,,. and 2d,,. states, is 
coupled to collective surface vibrations of the 
core. Good agreement is obtained with the known 
levels in these nuclei using a reasonable value 
for the coupling parameter. To obtain the agree- 
ment, the d,.-S,,. splitting, €, must be regarded 
as a function of neutron number. The manner in 
which € varies, as found in the intermediate- 
coupling calculation, is compared with the pre- 
dictions of the pairing correlation theory orig- 
inally introduced in connection with supercon- 
ductivity. Agreement as to the general trend is 
found. This may be regarded to some extent as 
an indication of the applicability of the pairing 
correlation theory to nuclear structure calcula- 
tions. 


PARITY -NONCONSERVING INTERNUCLEON 
POTENTIALS. R. J. Blin-Stoyle, Department 
of Physics and Laboratory of Nuclear Science, 
Massachusetts Institute of Technology, Cam- 


bridge, Massachusetts (Received January 22, 
1960). 


The general form that a parity-nonconserving 
internucleon potential must take because of in- 
variance requirements is obtained. A detailed 
calculation is then made of the parity-noncon- 
serving potential arising from a self-interacting 
current description of weak interactions. If the 
polar vector part of the current Ju") is con- 
served, then parity nonconservation of the order 


1 part in 10” ($~10~) is to be expected in nuclear 


Processes. Failure to observe such an effect 
would indicate either that J Wa is not conserved 





or that the self-interacting current description 
is incorrect. 


“On sabbatical leave from the Clarendon Laboratory, 
Oxford, England. 


SPINS AND DECAY MODES OF CERTAIN 
NEUTRON -DEFICIENT SILVER ISOTOPES. 

O. Ames, A. M. Bernstein, M. H. Brennan, 

R. A. Haberstroh, and D. R. Hamilton, Palmer 
Physical Laboratory, Princeton University, 
Princeton, New Jersey (Received January 22, 
1960). 


Isotopically enriched Pd foils were bombarded 
with protons to determine the origin of the 1.2- 
hr activity which was previously discovered in 
an even-A Ag isotope. The only appropriate 
activity was found in Ag'™ which was observed 
to have a (69+3)-min predominantly K -capture 
activity with intense y» rays of 550, 764, and 
920 kev. Using the atomic beam magnetic reso- 
nance method it was verified that the 69-min 
activity has J/=5. The /=2 resonance in Ag’™, 
when counted in an x-ray detector, was found 
to decay with a 69-min half-life with a small 
admixture of a 27-min component. When viewed 
in a 8 counter only the 27-min component is 
observed. These characteristics of the /=2 
resonance can be explained by placing the /=2 
level above the /=5 level with an appreciable 
amount of isomeric transition. This interpreta - 
tion is supported by the observation of feeding 
in the 920-kev y ray which only occurs in the 
decay of the /=5 state. 

Further work with isotopically-enriched Pd 
foils showed y rays of 120 and 150 kev, and, 
tentatively, 260 kev which have been assigned 
to the decay of 59-min Ag’®*. A positive identi- 
fication of a (15+ 2)-min activity in Ag’™ has 
been made. 


SOME CROSS SECTIONS FOR THE PRODUCTION 
OF RADIO-NUCLIDES IN THE BOMBARDMENT 
OF C, N, O, AND Fe BY MEDIUM ENERGY 
PROTONS. Masatake Honda and Devendra Lal,* 
Scripps Institution of Oceanography and School 

of Science and Engineering, University of Cali- 
fornia, La Jolla, California (Received November 
9, 1959; revised manuscript received March 14, 
1960). 


A number of nuclide formation cross sections 
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have been measured, using the Berkeley 184-in. 
cyclotron, to assist in the interpretation of the 
data on cosmic-ray-produced nuclides in the 
atmosphere and in iron meteorites. 

Cross sections of H*® and Be’ have been meas- 
ured in bombardments of organic targets con- 
taining nitrogen and oxygen by protons of energy 
225-730 Mev. Semicarbazide (CH,N,O) targets 
were used to obtain cross sections in air nuclei. 
The targets were prepared by mixing with a few 
percent of aluminum dust to permit reliable 
monitoring of the beam. Polyethylene, aluminum 
lactate, and Plexiglas targets provided elemen- 
tary cross sections in carbon and oxygen. 

The cross sections for the production of the 
long-lived isotopes C1°°(3 x 10° yr) and Al?*(8 x 10° 
yr) at 730 Mev, and of a number of short-lived 
radionuclides at 500 and 730 Mev, in iron bom- 
bardment by protons have been measured. These 
data and those of earlier workers suggest some 
modifications in empirical relations used for pre- 
dicting spallation cross sections in the case of 
nuclides close to stability. 


*on leave from the Tata Institute of Fundamental 
Research, Bombay, India. 


COSMIC-RAY PRODUCTION RATES OF Be’ IN 
OXYGEN, AND OF P*, P*%*, s** IN ARGON AT 
MOUNTAIN ALTITUDES. Devendra Lal, * 
James R. Arnold, and Masatake Honda, Scripps 
Institution of Oceanography and School of Science 
and Engineering, University of California, La 
Jolla, California (Received November 9, 1959; 
revised manuscript received March 14, 1960). 


The production rates of radio-isotopes P*, 
P**, and S** in argon, and of Be’ in oxygen have 
been measured by exposing argon and water to 
cosmic rays at mountain altitudes for periods of 
two to four months during 1959. The measured 
values at A=51°N, atmospheric depth 685 g cm~? 
are 7.6x10~°, 6.2x107*, 1.4x107° atom (gram 
argon)~ sec™ of P**, P**, and S*, respectively, 
and 9.0x10~° Be’ atom (gram oxygen)“ sec™. 

Isotope production rates for all regions in the 
atmosphere have recently been calculated by 
Lal et al. The measured production rates are 
higher than the calculated rates by factors of 
1.1, 1.8, 1.4, and 1.8 in the case of radio- 
isotopes Be’, P**, P**, and S**, respectively. In 
this comparison, account has been taken of the 
fact that cosmic-ray intensity has decreased by 
about 15% since 1948, the time period to which 
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the calculations apply. The measured production 
rates in oxygen are assumed to apply to air since 
cross sections for Be’ formation have been found 
to be the same in nitrogen and oxygen over most 
of the energy region of interest. 

The procedure used by Lal et al. yields a fairly 
accurate picture of the variation in production 
rates with altitude and latitude in the atmosphere. 
Isotope production rates in all regions of the at- 
mosphere can, therefore, be obtained by nor- 
malizing their calculations at the points where 
our measurements have been made. The cal- 
culated production rates of these isotopes in the 
troposphere, and in the stratosphere correspond- 
ing to observed cosmic-ray intensity during 1948. 
49, are given. The available data on the con- 
centration of the isotopes in rain-water and their 
averaged yearly deposition rates are compared 
with their revised production rates. From such 
a comparison, more definite conclusions can be 
drawn than hitherto possible in view of the more 
accurate knowledge of isotope production rates. 






















*on leave from the Tata Institute of Fundamental 
Research, Bombay, India. 





INELASTIC SCATTERING OF HIGH -ENERGY 
PROTONS EXCITING A COLLECTIVE LEVEL 
OF THE NUCLEUS. K. Nishimura, Department 
of Physics, Rutgers University, New Brunswick, 
New Jersey (Received October 9, 1959). 





ae ee lCU lC CU Cl eel 


An analysis is made of the angular distribution 
and polarization of 185-Mev protons inelastically 
scattered by carbon. A rough calculation using 
the distorted-wave Born approximation shows, 
even though the quantitative agreement is poor, 
that the 4.4-Mev level of C’* may be interpreted 
as a collective state. 


1-NUCLEON PHASE SHIFTS IN THE ENERGY 
RANGE 350 TO 600 Mev. W. D. Walker, J. Davis, @) 
and W. D. Shephard, University of Wisconsin, 
Madison, Wisconsin (Received January 27, 1960). By 


This paper contains the results of attempts to 
obtain a set of phase shifts valid in the energy 
region just above the 3-3 resonance. 

We have found it possible to extend the family 
of phase shifts found by the workers at Dubna 
at an energy of about 300 Mev. The features of 
this set of phase shifts are that a,, a@,, are posi- 
tive and a@,, is negative and small. We find in 
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dition that 5,, is positive and 6,, negative. At 
00 Mev the results are consistent with a reso- 
gace in the d,,. state. This resonant state, how- 
wer, decays a sizable fraction of the time into 
jfinal state with two 7 mesons. There are in- 
ications that at energies of 400 to 500 Mev 

nost of the single pion production comes from s 
ad py. States. 


{RAY YIELDS FROM ».-MESONIC ATOMS. 

y. A. Ruderman, Department of Physics, Uni- 
versity of California, Berkeley, California 
Received December 21, 1959). 


The interesting suggestion has been made that 
the rapid drop in the yield of mesonic K x rays 
inthe light elements may be associated with the 
apture of » mesons into the metastable 2s state. 
The mechanisms for making transitions from the 
isto the 1s state and from various p states into 
the 25 state have been investigated in detail for 
Li, Be, and B. It is found that the paradoxical 
reduction of K x rays remains unexplained. (1) 
Stark mixing of the mesonic 2s and 2p states by 
the electric fields of the atomic electrons allows 
“mixed” Auger -radiative transitions to the 1s 
state to compete favorably with radiationless 
transitions. These mixed transitions give a 
high-energy x ray and a relatively negligible 
(10-50 ev) electronic excitation and so contrib- 
wte to the observed K,, yield. (2) Even if the 
above “mixed” transitions are ignored, there 
isno mechanism which gets a large fraction of 
i mesons into the 2s state that at the same time 
does not violently contradict both theoretical 
estimates and observed K x-ray yields from light 
t‘mesonic atoms. 


MLARIZATION IN PION-NUCLEON SCATTER - 
NG AND THE SECOND AND THIRD PION- 
WCLEON RESONANCES. Michael J. Moravcsik, 
lawrence Radiation Laboratory, University of 
Glifornia, Livermore, California (Received 
hmary 25, 1960). 


The polarization of the recoil nucleon in pion- 
weleon scattering is studied from the point of 

liew of providing a means of distinguishing 

iong the various angular momentum assign- 
uents proposed for the higher pion-nucleon reso- 
lances. It is shown that polarization in this re- 
«ction is just as useful a guide as polarization 


in photoproduction. In particular, a measure- 
ment of the polarization at and in the neighbor - 
hood of 90° in the energy range between the first 
and third resonances should give a fairly con- 
vincing verification of one or the other of the 
prevailing assignments. 


EXAMPLE OF A SOLUBLE FIELD THEORY 
WITH FINITE CHARGE RENORMALIZATION. 
Herbert M. Fried, Department of Physics, Uni- 
versity of California, Los Angeles, California 
(Received January 11, 1960). 


A soluble field theory suggested by the Lee and 
Machida models is described in which coupling 
constant renormalization arises from a dressed 
boson and is finite if the contributing fermions 
are assumed nonrelativistic. For the unrenor- 
malized charge to be real, the renormalized 
charge must satisfy a certain inequality depend- 
ing on the boson and fermion mass ratios; if 
this inequality is violated a single boson ghost 
state occurs, as expected. 


FOLDY TRANSFORMATION IN THE PION- 
HYPERON SYSTEM. Akira Kanazawa, Depart- 
ment of Physics, Purdue University, Lafayette, 
Indiana (Received December 15, 1959). 


A unitary transformation, which plays the 
same role as the Foldy transformation in the 
pion-nucleon system, is constructed for the case 
where the pion interacts with both 2 and A hy- 
perons through y, couplings. The transformation 
function and the transformed Hamiltonian are 
very similar to those of the Foldy transforma- 
tion, in spite of the complexity of our system in 
isotopic spin space. The application to practical 
problems is not considered in this paper. 


GROUND-STATE ENERGY OF A MANY-FERMION 
SYSTEM. Il. J. M. Luttinger, University of 
Pennsylvania, Philadelphia, Pennsylvania, and 

J. C. Ward, Carnegie Institute of Technology, 
Pittsburgh, Pennsylvania (Received December 

7, 1959). 


The perturbation series for the ground-state 
energy of a many-fermion system is investigated 
to arbitrary order for the “isotropic” case. 
This is the case of over-all spherical symmetry, 
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both in the interaction and in the unperturbed 
single-particle energies. It is shown that for 
spin one -half fermions the Brueckner -Goldstone 
perturbation series is valid to all orders in the 
perturbation. For spins greater than one-half 

it is in general incorrect even in the isotropic 
case, unless the interactions are spin independ- 
ent. 

The discussion to arbitrary order in the inter - 
action is carried out by means of a Feynman- 
like propagator formalism, which is developed 
in detail. 


FORMAL PARADOX IN QUANTUM ELECTRO- 
DYNAMICS. Herbert M. Fried, University of 
California, Los Angeles, California (Received 
September 14, 1959; revised manuscript re- 
ceived November 30, 1959). 


The formal paradox concerning the vanishing 
of the photon self-mass, obtained by a mathe- 
matically improper manipulation, is examined 
in Kallén’s formulation of electrodynamics. It 
is suggested that this difficulty can be removed, 
and the formal manipulation retained, by a re- 
gularization of the Heisenberg operators. An 
alternate method of obtaining the spectral func - 
tion of the photon commutator is described, and 
a possible consequence of regularization, in 
connection with the proof of the renormalization 
constants’ divergence, is briefly discussed. 


GENERALIZATION OF QUANTUM MECHANICS. 
T. E. Phipps, Jr., Massachusetts Institute of 
Technology, Cambridge, Massachusetts (Re- 
ceived May 7, 1957; revised manuscript received 
January 20, 1960). 


The possibility of generalizing quantum mech- 
anics in such a way as to retain its predictive 
results, while comprehending additional solutions, 
is examined. It is found that this can be done 
through a perfected formal correspondence with 
Hamilton-Jacobi mechanics, by which one is led 
to consider generalizations of the Heisenberg 
postulate of the form Pp4j - Pp =S5.,, where S 
is a quantum analog of Hamilton’s principal 
function. The formalism is shown to be equiv- 
alent to a simple change in Hamiltonian, with 
transformed momentum operators satisfying 
conventional commutation relations, and with 
an additional relationship involving formal ana- 
logs of the classical “initial constants” adjoined. 
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A particular choice of S (=h/i) leads to a theory f # 
identical with wave mechanics apart from a 
constant (unobservable) phase factor on the ware | ! 
function. The fact that S may possess other, 
nonconstant values, demonstrated by a specific 
example, suggests the ability of the mechanical | | 
equations to describe a broader class of physica| 
states than has hitherto been investigated. | 


PION-PION SCATTERING IN THE ¢* THEORY, 
M. Baker and F. Zachariasen, Institute of Theo- 
retical Physics, Department of Physics, Stan- 

ford University, Stanford, California (Received 
December 21, 1959). 


Pion-pion scattering has been calculated using 
the determinantal method, assuming that a rela- 
tivistic (A/4)(9;9;)? coupling is responsible for 
the interaction. The scattering amplitude for the 
individual partial waves is expressed as a ratio 
of two power series and terms through A° have 
been kept in each series. Numerical results for 
the S and P waves have been obtained. is ad- 
justed by attempting to fit the electromagnetic 
structure of nucleons. The best value of A ob- 
tained by this fit is unfortunately so large that 
the validity of the determinantal approximation 
is doubtful. 


MOMENT OF INERTIA OF SUPERFLUID MANY- 
FERMION SYSTEMS. Ronald M. Rockmore, 
Brookhaven National Laboratory, Upton, New 
York (Received December 15, 1959). 


The effects of possible superfluidity on the 
cranking moment of a large many-fermion sys- 
tem moving under periodic boundary conditions 
are investigated within the framework of the 
theory of superconductivity recently formulated 
by Bogolyubov. The Hamiltonian is initially sub- 
jected to Bogolyubov’s general unitary trans- 
formation. The collective excitations of the 
fermions are then considered in the usual pair 
approximation; the appropriate cranking terms 
are linear in the boson pair operators. On per- 
forming a unitary transformation which trans- 
forms away these linear terms, one obtains an 
expression for the moment of inertia of the sys 
tem which includes both the effects of possible 
superfluidity and collective excitation. This ex 
pression, by virtue of its stationary property with 
respect to arbitrary variations in the amplitude 
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associated with the latter unitary transformation, 
is then utilized as a variational principle for the 
moment of inertia. For the normal state, the 
result previously obtained by the author, that the 
moment of inertia has the rigid value, is re- 
derived in more compact form. For the super - 
fluid state, one finds that collective excitations 
effect a marked increase in the superfluid mo- 
ment at intermediate coupling strengths although 
the resulting moment is still quite small com- 
pared to the rigid value. In contrast to the normal 
state case, where particle-hole pairs play a major 
role, this increase is almost entirely due to ex- 
citations consisting of particle-pairs or hole- 
pairs. The precise magnitude of the apparent 
resonance in the moment produced by the d-wave 
part of the cranking interaction is dependent to 
some extent on the features of the particle- 
particle potential which leads to the superfluid 
state. Variational expressions for the moment 
are exhibited for both Yukawa and delta-function 
shell potentials. These results are identical 

in charged and neutral Fermi systems. A cal- 
culation of the cranking moment at finite tem- 
peratures is presented in an Appendix along with 
an interpretation of it in terms of Bardeen’s two- 
fluid model of superconductivity. 


PUPPI-STANGHELLINI DISCREPANCY. H. P. 
Noyes, Lawrence Radiation Laboratory, Uni- 
versity of California, Livermore, California, 
and D. N. Edwards, Nuclear Physics Research 
Laboratory, University of Liverpool, Liverpool, 
England (Received December 21, 1959). 


In order to establish the extent of the dis- 
agreement between the pion-nucleon forward - 
scattering-amplitude dispersion relations and 
experiment in a statistical sense, the uncer- 
tainty in the dispersion integrals and S-wave 
scattering lengths is systematically included in 
the analysis. To accomplish this, we fit the total 
ross sections below 335 Mev by a Chew-Low 
P-wave resonance, phenomenologically modified, 
and calculate the error matrix for the param- 
eters. Our fit to the total cross sections is 
statistically at least as good as the Anderson 
darameterization used in previous work. Ig- 
toring forward scattering amplitudes above 220 
Mev because of D-wave uncertainties, we still 
find that there is less than 44% probability that 
the published data are compatible with a unique 
value for the pion-nucleon coupling constant f?, 


and that no adjustment of the S-wave scattering 
lengths can remove the discrepancy. However, if 
the 1” forward-scattering amplitudes measured 
by Ashkin et al. at 150 and 170 Mev are aban- 
doned in favor of the values recently obtained by 
Kruse and Arnold at 130 and 152 Mev, the prob- 
ability rises to 47.2% for our parameterization, 
or 8.6% for the Anderson parameterization. 

Cini et al. have pointed out that the conventional 
analysis of the low-energy data to obtain the S- 
wave scattering lengths does not satisfy cross- 
ing symmetry, and a reanalysis by Hamilton and 
Woolcock gives a, = -0.083, a_=0.088, rather 
than the conventional values of -0.110 and 0.077. 
We obtain some independent evidence in support 
of this conclusion by using the dispersion rela- 
tions to determine a,,a_, and f ? simultaneously. 
We find: 





Energy 
dependence ha a, a_ 
Anderson 0.075+0.018 -0.086 +0.025 0.071 +0.020 
Modified 


Chew-Low 0.086 +0.019 -0.101+0.026 0.085 +0.020 


It is clear that a better theoretical description 
of the energy dependence of the total cross sec- 
tions will be required before further progress 
can be made on this problem. 


LEVEL STRUCTURE OF NUCLEAR MATTER 
AND LIQUID He®. K. A. Brueckner and Toshio 
Soda, University of California, La Jolla, Cali- 
fornia, Philip W. Anderson, Bell Telephone 
Laboratories, Murray Hill, New Jersey, and 
Pierre Morel, French Embassy, New York, New 
York (Received January 15, 1960). 


Using the K matrix as computed in the study 
of nuclear matter and liquid He® as the effective 
interaction at the Fermi surface, the possible 
superfluidity of these systems has been investi- 
gated. A theory of the cooperative phenomenon 
valid for particle-particle interaction in states 
of arbitrary angular momentum has been devel- 
oped following the methods of Bardeen, Cooper, 
and Schrieffer. It is found for states of relative 
angular momentum other than /=0 that the parti- 
cle pairs must be correlated with respect to an 
arbitrary direction in the medium. As a result 
the change of structure of the Fermi surface 
is angularly dependent. An energy gap does not 
occur other than for /=0, the particle excitation 
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energy vanishing for certain orientations of the 
momentum. It is also shown that the specific 
heat shows a discontinuity at the transition tem- 
perature, but somewhat different from the case 
of 2=0. 

Application of these results to liquid He* shows 
that the cooperative effects arise from the in- 
teraction in the state with /=2, and that the 
transition temperature is at about 0.1°K. In 
nuclear matter the 'S, interaction is very weak 


and probably repulsive at the Fermi surface, 
and the attractive D, interaction gives a negli- 
gible energy shift. The *S, interaction is attrac. 
tive and in nuclear matter gives a few tenths of 
an Mev energy gap. These results suggest that 
in finite nuclei, with pairing of identical nucleons 
in the same shell, the cooperative effects are 
not strictly analogous to those in nuclear matter 
but instead are closely associated with the finite 


level spacing. 














